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PREFACE 


Although  there  is  a  significant  amount  of  reference  material  available  for  the  computer 
models  used  by  Air  Force  weather  forecasters,  there  is  no  single  reference  to  all  the  models. 
Ihis  techmcal  note  provides  that  reference  by  listing,  describing,  and  comparing  the  computer 
models  used  by  AFGWC  and  NMC. 

Because  of  the  complexity  of  atmospheric  processes,  niunerical  analysis  and  forecast  models 
use  approximations  based  on  theoretical  knowledge,  operational  considerations,  and 
available  computer  power.  The  strengths  and  weaknesses  of  the  varioxis  models  depend  on 
the  type  of  model  and  the  extent  of  the  approximations  that  are  used.  Weather  forecasters 
should  understand  these  strengths  and  weaknesses,  which  for  most  models  are  by  now  well 
documented. 

This  technical  note  provides  Air  Force  forecasters  and  others  with  basic  descriptions  of  the 
most  widely  used  analysis  and  forecast  models.  It  is  more  of  an  operational  guide  than  a 
comprehensive  models  textbook.  We  have  not  included  the  mathematical  equations  around 
which  the  models  are  built,  but  have  tried  to  provide  answers  to  such  questions  as,  "What 
is  a  model?",  "How  do  models  work?",  "What  are  the  strengths  and  weaknesses  of  the  models 
used?",  and  "Which  ones  apply  to  me?". 

Interested  readers  should  also  see  the  following:  AFGWCV87-001,  AFGWC  Cloud  Forecast 
Models',  AFGWC/88-001,  The  AFGWC  Automated  Real-Time  Cloud  Analysis  Model, 
AWS/86-001,  AJ^WC's  Advanced  Weather  Analysis  and  Prediction  System  (AWAPS)',  and 
AFGW<y79-004,  The  AFGWC  Automated  Analysis/ Forecast  Model  System. 

The  author  wishes  to  thank  the  following  for  their  help:  Major  Steven  Sycuro,  Capts  David 
Knapp  and  Michael  Fontaine,  and  Mr  Jay  Albrecht. 
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Chapter  1 

WHAT  iS  A  COMPUTER  MODEL? 


1.1  The  Model  Explained.  An 

atmospheric  computer  model  is  a  set  of 
mathematical  equations  that  describe  the 
movement  of  the  atmosphere  and  present  a 
graphic  view  of  current  and  forecast 
atmosphmc  conditions.  Ihe  first  successful 
nm  of  such  a  model  took  place  in  1950  at 
Princeton  University,  where  it  took  24  hours 
to  make  a  24*hour  forecast.  Model  forecasts 
are  not  perfect  for  the  following  reasons: 

•  We  have  an  incomplete  understanding  of 
the  atmosphere. 

•  Scientists  have  not  been  able  to  describe 
all  atmospheric  processes  mathematically. 

•  The  mathematical  equations  that 
describe  atmospheric  motions  do  not  provide 
exact  solutions;  they  must  be  approximated. 

•  Observations  are  geographically 
scattered;  in  some  areas,  they  are 
nonexistent.  Because  information  for 
data'sparse  areas  must  be  interpolated  from 
surrounding  gridpoints,  there  can  be  errors 
in  gridpoint  data. 

•  Computer  memory  size  and  central 
processing  unit  (CPU)  speed  are  limited. 
As  a  result,  only  a  limited  range  of 
atmospheric  features  can  be  predicted.  The 
smallest  feature  the  Air  Weather  Service 
Primitive  Equation  (AWSPE)  model  (no 
longer  in  use)  could  forecast  with  its  381-km 
grid  mesh  was  1,524  sq  km,  or  an  area 
about  the  size  of  the  Midwest. 


1.2  Improving  the  Model.  Computer 
models  can  be  improved,  but  to  do  so  is 
expensive  and  time-consuming.  For 
example,  while  forecast  accurate  can  be 
improved  10-15%  by  halving  the  grid  space, 
it  would  take  four  times  as  long  to  run  such 
a  model  on  the  same  computer.  The  time 
step  (the  time  it  takes  for  the  fastest  model 
wave  to  move  firom  one  gridpoint  to  another) 
would  have  to  be  halved,  increasing  the 
CPU  requirement  by  a  factor  of  two.  Ihe 
combined  resxilt  would  be  an  increase  by  a 
factor  of  eight;  that  is,  if  the  AFGWC  Global 
Spectral  Model  (GSM)  took  20  minutes  to 
run  normally,  it  would  take  160  minutes  to 
run  at  two  tunes  the  resolution  on  the  same 
computer.  This  alone  makes  increasing  the 
grid  resolution  by  two  an  ineffective  way  to 
increase  forecast  skill. 

1 .3  Problems.  The  major  problems  with 
numerical  weather  prediction  are: 

1.3.1  Physics  (about  30%).  Physics  problems 
stem  fix>m  our  lack  of  a  complete 
tmderstanding  of  the  atmosphere  and  the 
resultant  need  to  approximate  some 
processes  because  computer  size  and  speed 
are  limited.  For  example,  turbulence  and 
friction  are  parameterized,  meaning  that  the 
calculation  of  subgrid  scale  processes  (those 
processes  that  are  too  small  to  be  detected 
on  the  grid)  are  approximated  on  a  grid- 
scede  basis. 

1.3.2  Initial  conditions  (about  25%)  are  the 
most  difficult  to  correct.  They  are 
problematic  because: 


*  Weather  observmg  statiom  are  located 
at  random,  rather  than  at  the  exact 
gridpoiota  of  the  model.  Observatioiis  are 
concentrated  in  some  aieas,  such  as  the 
United  States  and  Btirope,  and  sparse  in 
othera,  such  as  Asia. 

a  Weather  observations  often  contain 
errors. 

•  Weather  observations  are  often 
inconsistent,  partly  because  of  the  methods 
used  to  collect  the  data  and  partly  because 
of  calculation  errota  (for  example,  RAOB 
and  satellite  temperatures  for  the  same 
point  may  be  different). 

1.3.3  Raaohition  (about  40%).  In  a  grid 
model,  the  smallest  wave  that  can  be  seen 
or  resolved  horizontally  is  two  grid  lengths-, 
that  is,  if  the  grid  space  were  200  km,  the 
smallest  resolvable  wave  would  be  400  km 
long. 

Aliasing.  Ihuring  the  forecast  process,  false 
waves  that  are  fotir  or  lesa  grid  lengths  long 
are  created  by  the  mathematical  terms 
within  the  model.  These  false  waves,  if  not 
removed,  will  cause  the  model  to  become 
imstable. 


Tima  Sb^.  The  time  step  in  a  model  must 
be  smaller  than  a  certain  value  or  the  model 
will  become  unstable  and  produce 
meaningless  restdis. 

PhaaaSpaad.  Gridpoints  in  most  numerical 
models  are  fixed,  both  horizontally  smd 
vertically.  Because  of  this,  small 
perturbations  occurring  during  the  forecast 
process  tend  to  move  more  slowly  than  they 
would  in  reality. 

Spaetral  v$.  Grfdpolnt  Modala. 

Gridpoint  mo(Ma  use  U-V  wind  components  to 
represent  the  vector  value  of  the  wind  at 
each  individual  gridpoint.  Wind  is  a  vector 
that  includes  both  speed  and  direction 
Wind  vectors  can  be  resolved  into  eastward 
(U)  and  northward  (V)  components. 
Analysis  and  forecast  models  resolve  wind 
vectors  into  west-to-east  (+U)  and 
south-to-north  (+V)  components. 

Spectra/  mocMs  use  sine  and  cosine  functions 
to  represent  weather  elements  such  as  wind 
and  temperature.  These  sine  and  cosine 
functions  go  tiurough  changes  in  amplitude 
during  the  forecskst  process.  Because  of  this 
special  characteristic,  spectral  models  do  not 
have  the  aliasing  and  phase  speed  problems 
mentioned  above. 


2 


Ch^er  2 
GRIDS 


2.1  AFGWC  Grids.  For  our  purposes,  the 
only  grid  ^jrpes  discussed  (unless  otherwise 
spewed)  are  those  in  use  at  AFGWC. 

2J2  Grid#  Ba«jc  to  Modots.  The  starting 
point  for  any  gridpoint-based  model  is  the 
grid.  The  atnu^phere  is  relatively  shallow 
in  die  vertical,  but  broad  horizontally.  This 
is  why  die  number  of  gridpoints  in  the 
horizontal  is  usually  greats  than  the 
number  of  gridpoints  in  the  vertical.  In  die 
horizontal,  the  gridpoints  are  equally 
spaced.  In  die  vertical,  points  are  "stacked” 
to  represent  those  areas  of  the  atmosphere 
that  are  commonly  analyzed  or  that  are  of 
significance  to  the  model  (e.g.,  860  mb,  7(X) 
mb,  500  mb).  Grid  spadngs  and  vertical 
stacking  are  chosen  in  such  a  way  as  to 
simplify  the  calculations  of  the  forecast 
model  and  to  provide  outputs  tailored  for  a 
specific  customer's  use. 

2.3  The  Grid  Described.  A  grid  is  a 
group  of  regularly  spaced  points  that 
represent  the  intersection  of  regularly 
spaced  parallel  and  perpendicular  lines'-see 
Figures  1  and  2.  The  standard  and  lowest 
resolution  grid  used  at  AFGWC  is  the  whole 
mesh.  Higher-resolution  grids  (or 
finer-mesh  grids)  use  whole-mesh  gridpoints 


as  a  starting  point  and  have  more  gridpoints 
in  the  horizontal.  For  example,  the 
half-mesh  grid  is  simply  a  grid  that 
contains  the  whole  mesh  gridpoints  and  one 
gridpoint  between  every  whole  mesh 
gridpoint.  Figure  3  is  an  example  of  the 
relationship  between  whole-mesh, 
half-mesh,  and  quartor-mesh  gridpoints. 


Figure  1.  An  Example  of  a  Grid.  A  grid  is 
simply  an  array  of  points  located  at  the 
intersections  of  uniformly  spaced  parallel 
and  perpendicular  lines  (from  Hoke  et  al., 
1986). 
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Figure  2.  An  Exampte  of  a  Grid 
Suparimpoaad  on  a  Map.  (From  Hoke  et  al., 
1^). 


2.4  Grid  Resolution.  Atmospheric 
motions  on  a  scale  smaller  dian  the  grid  in 
use  cannot  be  represented  conrectly  because 
wavelengths  must  be  equal  to  or  greater 
than  twice  the  distance  between  two 
adjacent  gridpoints.  Also,  a  wave  is 
genmnlly  handled  poorly  if  its  wavelength  is 
less  than  four  gridpoints  wide.  Since  many 
atmospheric  events  deal  with  motions 
smaller  than  1,6(X)  km-*the  distance  covered 
by  four  successive  whole-mesh 
grii^ints'-grids  with  finer  resolution  were 
developed.  These  include  the  half-mesh, 
quarter-mesh,  eighth-mesh,  and 
sixty-fourth-mesh  grids.  As  their  names 
imply,  the  distance  between  gridpoints  on 
ihe  finer  meshes  is  defined  relative  to  the 
whole-mesh. 
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Figura  3.  Tha  Raiationship  Bataraan 
Wiiaia<Maah,  Haif*Maah,  and  Quartar*Maah 
Grtdpmnta,  Note  that  each  wbole-mesh 
point  is  a  half-mesh  and  quarter-mesh  point, 
as  well.  Similarly,  each  half-mesh  point  is 
also  a  quarter-mesh  point  (from  Hoke  et  al., 
1985). 

2.5  The  Whole-Mesh  Grid.  Whole-mesh 
gridpoints  are  exactly  381  km  apart  at 
60®  N  relative  to  the  Earth's  surface.  This 
means  that,  because  of  the  Earth's 
curvature,  resolution  decreases  toward  the 
poles  and  increases  toward  the  Equator.  At 
half-mesh,  the  gridpoints  are  190.5  km 
apart.  Higher  resolution  grids  have 
successively  smaller  distances  between 
gridpoints.  An  example  of  the  whole-mesh 
grid  is  shown  in  FigU3:e  4. 
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Figur*  4.  The  Northern  Hemisphere  Whole>Mesh  Grid.  The  grid  spacing  is  displayed  in  the 
upper  left-hand  comer.  Note  that  the  grid  extends  beyond  the  boundaries  of  this  map.  The 
1^  indexing  convention  is  also  displayed.  This  convention  allows  for  labeling  of  every 
gridpoint.  I,  and  J,  coordinates  represent  the  whole-mesh  grid  scale.  I2  and  J2  represent  the 
half-mesh  grid,  and  so  on  (firom  Hoke,  et  ai.,  1985). 
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2.6  Hemispheric  Data  Displays.  The 

following  are  two  examples  of  gridded 
hemispheric  data  displays  (dn*  and  paper 
facsimiles): 

•  Hie  first  uses  all  whole-mesh'beised 
gridded  data  firom  the  poles  to  the  Equator. 
Gridpoints  beyond  the  Equator  are  ignored. 

•  The  second  is  the  octagon,  shown  in 
Figure  5.  An  example  of  gridded  data  in  an 
octagonal  display  is  the  Cloud-Free 
information  on  AFGWC  Horizontal  Weather 
Depiction  (HWD)  products.  Data  in  this 
display  extends  from  the  pole  to  the  borders 
of  the  octagon.  All  data  beyond  the 
octagon's  borders  is  ignored.  Limitations  of 
gridded  data  displays  are  generally  based  on 
computer  processing  capabilities  and  the 
printing  equipment  available  at  the  time  the 
model  was  designed. 

2.7  Stationary  Window  Grids  use 

windows,  or  subsets  of  the  whole-mesh  grid, 
to  cover  specific  lureas  of  the  globe  such  as 
the  North  American,  European,  and  Asian 
continents. 

2.8  The  Satellite  Global  Database 
(SGDB)  Grid.  This  64th-mesh  grid 
contains  16,777,216  points  in  the  Northern 
Hemisphere.  It  is  used  in  the  SGDB  display 
on  the  Satellite  Data  Handling  System 
(SDKS),  with  a  resolution  of  about  5  km. 


2.9  Grid  Types.  There  are  three  grid 
types  in  use  at  AFGWC:  polar 

stereographic,  Mercator,  and  latitude- 
longitude. 

Polarstereogrephic  (PST)  Gride.  The 
Northern  (or  Southern)  Hemispheric 
Whole-mesh  Reference  Grid  (sometimes 
called  the  "Northern  Hemispheric 
Whole-Mesh  Super  Grid")  is  a  66  X  65  point 
whole-mesh  grid  based  on  the  PST  map— see 
Figure  6.  The  domain  for  all  PST  gri^  are 
centered  on  the  hemisphere  (Northern  or 
Southern)  and  include  the  entire 
hemisphere.  The  borders  of  the  grids 
partially  extend  into  the  opposite 
hemisphere,  but  data  in  these  overlap  areas 
are  usually  ignored.  PST  grids  are  available 
in  the  half-,  quarter-,  eighth,  and  64th- 
mesh  resolutions.  They  are  used 
extensively  at  AFGWC  for  various  weather 
support  functions. 

Mercator  Grida.  These  grids  cover  an  area 
from  50®  N  to  50®  S  for  tropical 
meteorological  needs.  They  are  based  on  the 
mercator  map  projection. 

•  The  Conventional  Tropical  Grid  is  used  for 
conventional  meteorological  elements  such 
as  temperature,  height,  and  wind-see 
Figure  6. 

•  The  Satellite  Global  Database  Tropical  Grid  is 
used  for  procssing  satellite  imagery;  see 
Figure  7. 
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Figure  5.  An  Example  of  a  Polar-Stereographic  (PST)  Grid  (Whole-Mesh)  with  the  Octagon 
Grid  Boundaries  Centered  on  the  Northern  Hemisphere.  The  dashed  lines  indicate  that  the 
data  on  the  octagon  grid  is  actually  stored  in  a  rectangular  database.  All  data  in  the  comers 
of  the  rectangle  (outside  the  octagon)  are  ignored.  The  reference  grid  for  the  octagon  is  the 
outer  whole-mesh  grid  (from  Hoke  et  al.,  1986). 
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Conv«ntionaf  Tropica!  Grid  is  an  Example  of  a  Mercator  Grid. 


T 


in 
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Orki*.  The  Global  atAFGWC;  an  example  is  shown  in  Figwe 
Applications  Database  (GADB)  grid  and  the  8.  Grid  spacing  on  both  is  2.5“  x  2.5“.  Both 
High  Resolution  Analysis  System  (HIRAS)  contain  10,585  points, 
grid  are  the  latitude-longitude  grids  in  use 
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Chapter  3 

THE  AFQWC  PRODUCTION  CYCLE 


3.1  The  Advanced  Weather  Analysis 
and  Prediction  System  (AWAPS)  was 

installed  at  AFGWC  in  1985,  and  formally 
implemented  in  June  1986.  It  produces 
nearly  all  of  AFGWC's  automated  forecasts 
except  for  cloud^\  AWAPS  comprises  four 
models:  the  Global  Spectral  Model  (GSM), 
the  First-Guess  model,  the  High  Resolution 


Analysis  System  (HIRAS)  model,  and  the 
Relocatable  Window  Model  (RWM).  Figure 
9  shows  how  the  GSM  smd  HIRAS 
processing  cycles  are  linked  to  form  part  of 
the  AFGWC  production  cycle,  a  series  of 
forecasts,  analyses,  and  reanalyses  that  runs 
at  6-hour  intervals. 


VALID  TIML  |Z| 


Figure  9.  The  HIRAS/GSM  Relationship  in  the  AWAPS  Production  Cycle  (from  Hoke  et  ai., 
198S). 
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3.2  The  AWAPS  Cycle.  Both  the  GSM 
and  First'Guess  cure  updated  every  6  hours; 
although  slightly  different,  they  provide 
similar  product  support.  The  HIRAS  uses 
the  First'Guess  data  fields  to  produce  the 
best  possible  initial  analysis  for  the  next 
6*hour  cycle.  The  OOOOZ  cycle  is  shown 
below  as  an  example: 

At  0100Z,  the  Fii^t'Guess  is  ready  for 
quality  control  (QC)  by  forecasters.  "Bogus" 
data  (this  will  be  explained  in  a  later 
chapter)  can  be  input  to  alter  the 
First'Guess  data  ffelds  and  improve  the 
analysis. 

At  0130Z,  the  HIRAS  surface  analysis  is 
run.  It  takes  10  minutes  to  lom  the  surface 
analysis  on  the  CRAY/X-MP  mainframe 
computer.  Satellite  soundings  are  anchored 
periodically  using  this  or  one  of  the  next 
three  sxurface  analyses.  Satellite  soimdings 
are  built  by  adding  thickness  values  to  the 
surface  HIRA.S  1,000  mb  analysis  ffeld. 

At  0230Z,  the  HIRAS  surface  analysis  is  run 
for  the  second  time.  It  is  also  used  to 
anchor  the  satellite  soundings  for  the  first 
upper-air  analysis. 

At  0300Z,  the  HIRAS  upper-air  analysis  is 
run  on  the  CRAY/3f/MP.  It  takes  20 
minutes. 

At  0400Z,  the  GSM  is  run  out  to  96  hours 
using  the  HIRAS  first  upper-air  analysis 


for  its  initial  conditions.  It  takes  1.5  hours 
to  run  the  model  on  the  CRAY/X-MP. 

At  0530Z,  HIRAS  does  its  third  surface 
analysis;  it  is  again  used  to  anchor  satellite 
sounding  data. 

At  0600Z,  a  second  upper-air  analysis  is  run. 
OOOOZ  data  will  continue  to  come  in  up  to 
this  point;  all  the  new  data  available  is 
essential  for  making  the  best  possible  OOOOZ 
analysis. 

At  0630Z,  the  First-Guess  model  makes  its 
6-hour  forecast  using  the  second  HIRAS 
OOOOZ  upper-air  analysis  for  its  inital 
conditions.  It  takes  30  minutes  to  run  out 
to  12  hours  on  the  CRAY/X-MP.  When 
finished,  the  cycle  begins  again. 

3.3  First-Guess  Accuracy.  In  order  to 
have  a  more  accurate  First-Guess,  the  model 
produces  a  forecast  every  6  hours.  This 
forecast  is  then  used  by  HIRAS.  A  6-hour 
cycle  provides  a  higher  quality  analysis  than 
would  a  12-hour  cycle;  the  analysis  is 
updated  more  often  euid  used  for  initializing 
the  next  forecast  run,  which  in  turn 
produces  a  better  forecast. 

3.4  The  Relocatable  Window  Model 
(RWM)  is  being  developed  for  inclusion  in 
the  AWAPS  production  cycle.  It  is  designed 
to  make  a  high-quality  and  high-resolution 
forecast  anywhere  in  the  world  at  any  time. 
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Chapter  4 

THE  HIGH  RESOLUTION  ANALYSIS  SYSTEM  (HIRAS) 


4.1  The  HIRAS  Model.  HIRAS  is  an 
analysis  model  used  primarily  to  initialize 
the  forecast  models.  It  is  extremely 
important.  HIRAS  deficiencies  can 
dramatically  alter  the  outcome  of  the 
forecast  models. 

4.2  HIRAS  Input  Data.  HIRAS  uses 
observations  from  land  stations,  ships, 
buoys,  aircsraft,  radiosondes,  and  satellites 
(see  Figures  10'15  for  typical  global 
distributions  of  the  various  observation 
types).  Although  HIRAS  applies  most  of  the 
data  available  to  it  to  its  internal  grid,  it 
cant  use  all  the  data  because  of  software 
and  hardware  limitations.  The  output  of  the 
First-Guess  model  is  essential  to  the  run  of 
the  HIRAS  model.  HIRAS  blends  its  own 
data  with  the  output  from  the  First-Guess 
model  (a  6>liour  forecast  valid  at  the  same 
time)  to  create  the  best  possible  analysis. 

4.3  Optimum  intcsrpolation  (Ol).  HIRAS 

does  its  analysis  using  a  technique  called 
"Optimiim  Interpolation,"  an  analysis 
technique  used  by  most  contemporary 
analysis  models.  Our  discussion  of  OI  is 
based  on  the  version  used  by  HIRAS. 
Basically,  there  are  three  things  OI 
considers: 

•  The  number  of  observations  available 
and  the  distance  between  them  end  the 
analysis  grid  point. 

•  The  precision  of  the  equipment  used  to 
take  the  observations. 

•  The  statistical  accuracy  of  the 
First-Guess  fields. 


4.3.1  Observation/Gridpoint  Distance.  The  use 
of  the  mean  distance  between  observations 
is  not  unique  to  01;  it  is  found  in  nearly 
every  other  numerical  analysis  scheme.  In 
very  simple  terms,  the  model  assigns 
weights  to  observations  surrounding  each 
gridpoint.  Each  observation  is  allowed  to 
influence  surrounding  gridpoints  by  the 
weight  it  carries;  weights  decrease 
exponentially  with  distance.  Observations 
in  the  vertical  are  also  assigned  weights. 
All  the  observations  and  their  weights  are 
applied  mathematically  to  the  First-Guess 
data  at  all  gridpoints  and  atmospheric 
levels  to  modify  the  First-Guess  data  as 
needed. 

4.3.2  Accuracy  of  the  Observing  Instrument 

The  model  is  able  to  distinguish  between  the 
various  observing  insii’uments  used.  Each 
instrument  type  is  assigned  a  statistically 
determined  expected  error.  The  lower  the 
expected  error,  the  more  weight  the 
observation  carries.  For  example,  a  600-mb 
temperature  from  a  radiosonde,  which  is 
probably  more  accurate  than  the  500-mb 
temperature  from  a  satellite  sounding,  is 
assigned  a  higher  weight. 


13 


Figure  10.  Typical  Distribution  of  Surface  Marine  Observations  at  1200Z  (from  Day,  1989). 


Figure  11.  Typical  Distribution  of  Aircraft  and  ASDAR  (Aircraft  to  Satellite  Data  Relay) 
Observations  at  1200Z  (from  Dey,  1989). 


Figure  13.  Typical  Distribution  of  Satellite  Temperature  Profile  Soundings  at  1200Z  (from 
Day,  1989). 


1*.  Tirpici  Oftlribulion  of  Cloud-T™ck«l  Wind  Obonvotlon.  of  1200Z  (from  Dey, 


250-mb  HiRAS  Error  Fields  as  Applied  to  the  First-Guess 
Stobie,  1986). 


4.3.3  First-Guess  Accuracy.  HIRAS 
produces  a  set  of  Reids  (called  "error  fields") 
that  indicate  how  accurate  the  analyses  a3% 
at  all  gridpoints.  Essentially,  the  more 
observations  there  are  to  influence  given 
gridpoints,  the  more  accurate  the  resulting 
analysis  will  be.  Figure  16  (opposite)  is  an 
example  of  a  250-mb  height  error  Reid. 
Note  that  error  values  are  low  ever  data-rich 
areas  such  as  North  America  but  high  over 
data-sparse  areas  such  as  the  Indian  Ocean. 
In  addition  to  indicating  analysis  accuracy, 
error  Reids  also  indicate  the  accuracy  of  the 
Furst-Guess. 

To  see  how  this  works,  let's  consider  two 
areas,  North  America  and  the  South  PaciRc. 
Since  there  are  a  large  number  of 
radiosonde  stations  there,  the  analysis  over 
North  America  is  considered  accurate.  But 
there  are  very  few  upper-air  stations  in  the 
South  PaciRc;  although  a  satellite  makes  an 
occasional  pass,  satellite  data  is  not  nearly 
as  reliable  as  radiosonde  data.  As  a  result, 
initial  conditions  in  the  South  PaciRc  are 
usually  poorly  handled  by  the  First-Guess 
model  and  the  First-Guess  over  the  South 
PaciRc  is  not  usually  as  accurate  as  over 
North  America. 

Since  there  is  this  difference  in  First-Guess 
accuracy  from  place  to  place,  HIRAS  01 
gives  more  weight  to  the  First-Guess  data 
when  compared  to  any  individual 
observation  in  data-rich  areas  like  North 
America  and  more  weight  to  individual 
observations  in  data-sparse  areas  like  the 
South  Pacific.  The  aggregate  of 
observations  in  data-rich  areas,  however, 
usually  esury  more  weight  when  applied  to 
an  individual  First-Guess  gridpoint.  The 
obvious  question  is  "Why  doesn't  the  HIRAS 
always  place  more  emphasis  on  observations 
from  data-rich  areas?" 


01  uses  a  combination  of  observations  and 
the  First-Guess.  The  better  the  First-Guess 
(and  it  usually  is  better  in  data-rich  areas) 
the  less  the  weight  given  to  individual 
observations.  Observations  also  contain 
errors.  A  radiosonde,  for  example,  might 
have  a  3®  C  temperature  error.  If  the 
First-Guess  had  an  expected  error  of  only  1° 
C,  it  wouldn't  make  sense  to  use  oiUy  the 
RAOB.  All  the  data  available  around  any 
given  gridpoint  is  used  unless  the  QC 
function  determines  the  data  to  be  bad.  01 
determines  how  much  each  piece  of  data  will 
affect  the  final  output  for  a  given  gridpoint. 

Although  it  doesn't  happen  often,  HIRAS 
may  not  recognize  dramatic  changes  in 
atmospheric  conditions  quickly.  The 
First-Guess  model  has  considerable  weight 
in  data-rich  areas.  If  the  observations  show 
a  dramatic  difference  and  the  First-Guess  is 
really  wrong,  only  a  portion  of  the  difference 
will  be  used  to  affect  the  outcome  of  the 
analysis.  In  rare  circumstances,  it  may  take 
a  couple  of  forecast  cycles  for  HIRAS  to 
recognize  rapid  changes. 

4.4  The  Two  HIRAS  Models.  HIRAS  is 

really  two  different  models:  upper-air  and 
surface. 

4.4.1  The  Upper-Air  Analysis  extends  from 
1,000  to  10  mb  and  is  used  to  initialize  the 
GSM.  HIRAS  directly  analyzes  only  five 
different  Reids:  geopotential  heights, 
U-component  wind  (West-East), 
V-component  wind  (South-North), 
temperature,  and  relative  humidity.  From 
these  five,  HIRAS  derives  thirteen  other 
Reids.  A  number  of  other  Reids  not 
produced  by  HIRAS  (e.g.,  thickness  and 
vertical  velocity)  are  produced  on  other 
computer  processors  using  the  HIRAS 
Reids.  Table  1  provides  a  list  of  all  the 
Reids  produced  by  HIRAS  or  derived  by 
HIRAS  outputs. 
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TABLE  1 .  U»t  of  floldt  analyzod  by  or  darivad  from  tha  HtRAS  upper-air  model.  Niunerous 
othw  fields  are  processed  through  algorithms  on  SDHS.  The  list  includes  the  levelVlayez« 
in  millibars  for  which  the  data  are  available. 


FflELO 

<lir«ciiy  awlyzcd 

Height/Layer  in  Millibars 

Itebrtiyc  Humidity 

Sfe,  1000,  eso,  700,  500, 400,  JOO 

Hc^ihts 

1000,  OSO,  700,  SOC,  400,  300,  3S0, 200,  ISO,  100,  70,  SO,  30,  20,  10 

Vsmpemture 

Sfo,  1000,  ISO,  700,  500,  400,  300,  2S0,  200,  ISO,  100,  70,  50,  30,  20, 10 

U-omipontiit  wbids 

ate,  1000,  850.  700,  S08,  400,  300,  250,  200,  ISO,  100,  70,  50,  30,  20,  10 

V-Compon«nt  wbid 

sfC,  1000,  050,  700.  500,  400,  300,  250,  200,  ISO,  100,  70,  58,  30,  20, 10 

FIELDS 

Helght/Uyer  in  MilUbart 

derived 

Tempeiature  Advection 

SO 

Streamline  Functions 

1000,  ISO,  700,  SOO,  400,  300,  250,  200, 150, 100,  70,  SO,  30,  20,  10 

spedfle  Humidity 

Sfo,  1000,  ISO,  700,  SOO,  400,  300 

Tropopaute  Haight 

• 

TropopauK  Pressure 

* 

Tropopause  Temperature 

e 

Sea  Level  Pressure 

Dewpoint 

Sfo,  1000,  ISO,  700,  500,  400,  300 

1  Precipitahic  water 

Sfo-OSO,  0 S0'70«,  700-500,  500-400,  400-300,  300-100 

1  vortkity  Advection 

•SO,  SOO 

I  vorticlty 

150,  500,  300 

1  "Packed  U-v  Winds 

SfC,  1000,  ISO,  700,  SOO,  400,  300,  2S0,  200, 150, 100,  70,  SO,  30,  20,  10 

1  ‘"Hours  Since  land 

sfc  (boundary  layer) 

*  Tropopause  data  varies  with  rest 

met  to  its  exact  position. 

1  ‘  *  Packed  U-V  Winds  •  the  individual  U  and  V  wind  components  contain  too  many  bits  of  information  and  must  be 

1  OMnpacted  for  SONS  to  use  them. 

1  "  *  Describes  the  amount  of  time  a  boundaiy  layer  parcel,  currently  over  water,  has  spent  over  water. 
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The  upper*air  analysis,  which  is  displayed 
using  a  2.5**  x  2.5°  latitude-longitude  grid,  is 
actually  run  on  the  Kurihara  grid  (shown  in 
Figtire  17)  to  save  computer  time.  The 
Kurihara  grid,  with  a  coaraer  (or  lower) 


resolution,  has  about  a  third  as  many 
gridpoints  as  the  display  grid.  The  upper- 
air  analysis  uses  a  spectral  transformation 
along  writh  01. 
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The  first  step  m  spectral  tramformation  is 
to  take  the  First-Guess  30-wave  spectral 
coeflhcients,  traxksform  them  into  gridpoint 
values,  and  place  them  on  the  coarser 
HIRAS  grid.  When  this  is  completed,  the 
model  performs  its  OI.  When  the  01  is 
completed,  the  model  transforms  its 
gridpoint  oaalysis  into  40-wave  spectral 
coefficients.  Although  the  analysis  is  done 
using  gridpoints,  the  conversion  to  spectral 
coefficients  allows  each  gridpoint  to  affect 
the  entire  analysis.  In  other  words,  a 
gridpoint  temperature  in  Ohio  could 
influent  a  gridpoint  temperature  in  the 
Caspian  Sea.  This  transformation  also 
allows  for  the  final  transformation  back  to 
the  2.5°  X  2.6°  grid.  Even  though  the 
analysis  appears  in  the  database  on  a  2.5° 
X  2.5°  grid,  the  true  analysis  resolution  is 
that  of  the  coarser  Kurihara  grid  that  was 
shown  in  Figure  17. 

Because  of  deficiencies  in  the  First-Guess 
model  above  100  mb,  HIRAS  uses  the 
previous  HIRAS  analysis  (6  hours  old)  for 
its  first-guess  instead.  The  model  forces  the 
winds  above  100  mb  to  be  geostrophic  and 
uses  only  horizontal  OI  above  100  mb"that 
is,  it  doesn't  look  up  or  down  to  interpolate 
data  at  higher  or  lower  levels.  As  a  result, 
the  analysis  above  100  mb  is  two-,  rather 
than  three-dimensional. 

4.4.2  The  HIRAS  Surface  Analyeie  derives 
1,000  mb  temperature  and  geopotential 
height  fields  from  the  surface  temperature 
and  pressure  fields.  It  is  used  primarily  to 
anchor  satellite  sounding  data,  which 
includes  a  great  deal  of  information  for  the 
different  mandatory  pressvire  levels,  but 
lacks  the  geopotential  height  of  each 
pressure  level;  therefore,  the  1,000-mb 
height  Belds  from  the  HIRAS  surface 
analysis  are  used  to  determine  the  height  of 
the  1,000-mb  pressure  level  in  the  satellite 
soundings.  The  rest  of  the  pressure-level 


geopotential  heights  can  be  determined  from 
this  value. 

The  surface  analysis  is  also  an  OI  analysis. 
It  differs  from  its  upper-air  counterpart 
because  it  performs  its  analysis  directly  on 
the  2.5°  X  2.5°  grid.  It  does  not  use  spectral 
transformation  or  internal  grids.  It  uses  the 
OI  error  fields  to  ensure  grid  data  accuracy, 
but  observations  are  weighted  more  heavily 
in  the  surface  analysis  than  in  the  upper-air 
analysis  throughout  the  hemisphere. 

4.5  Quality  Control.  Since  observation 
errors  are  routine,  manual  and  automated 
quality  control  (QC)  must  be  performed. 
Bad  observations  must  be  corrected  or 
discau^ed  in  order  to  prevent  the  model  from 
providing  a  bad  analysis.  QC  is 
accomplished  in  two  ways: 

4.5.1  Manual  QC  is  performed  by  AFGWC 
forecasters,  who  delete,  change,  or  create 
observations  based  on  known  information. 
This  forecaster-edited  data  is  known  as 
"bogus"  data. 

4.5.2  Automated  QC.  An  automated  QC 
system  is  built  into  the  model.  The  system 
can  recognize  a  bad  observation  (e.g.,  a  990- 
mb  sea-level  pressure  within  a  1,046-mb 
high)  and  discard  that  observation.  There 
are  two  steps  to  this  process: 

•  Gross  CheckB.  Observations  are 
compared  to  adjacent  gridpoint  data  from 
the  First-Guess  model  output.  If  there  is  a 
large  enough  disagreement,  the  observation 
is  discarded.  If  it  just  barely  passes  the 
gross  check,  it  is  held  over  for  the  second,  or 
"buddy"  check. 

•  Buddy  Check.  An  observation  that  barely 
passes  the  gross  check  is  compared  with 
stations  nearby;  if  the  observation  disagrees 
significantly,  it  is  thrown  out. 
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Chapter  5 

THE  GLOBAL  SPECTRAL  MODEL  (GSM) 


5.1  Introduction.  The  GSM  (AFGWC's 
primary  forecast  model)  was  designed 
primarily  as  an  aviation  model.  There  are 
two  GSM  models  in  use  at  AFGWC;  the 
GSM  proper  and  the  First-Guess.  The  GSM 
and  l^e  First  Guess  are  run  every  6  hours. 
The  GSM  proper  produces  forecasts  out  to 
96  hours  and  First-Guess  out  to  12  hours. 
Although  the  GSM  provides  a  wide  variety 
of  forecast  fields  (see  Tables  2  and  3),  its 
main  purpose  is  to  provide  height,  wind,  and 
temperature  forecasts  for  aviation 
customers.  The  First-Guess  model  differs 
slightily  fiom  the  GSM.  The  material  in  this 
chapter  includes  the  First-Guess  forecast 
model  unless  otherwise  specified. 

5.2  Spectral  Forecasting.  A  true 
spectral  representation  iises  mathematical 
combinations  of  sines  and  cosines  to 
represent  patterns  of  meteorological  fields. 
That  is,  a  spectral  equation  contains  a  series 
of  terms,  each  having  either  a  sine  or  cosine 
term  as  a  factor.  Becaiise  sines  and  cosines 
look  like  waves,  spectral  representations  are 
often  discussed  in  those  terms;  the  GSM, 
for  example,  is  a  40-wave  model,  even 
though  it  uses  many  more  than  40  sine  and 
cosine  terms.  In  fact,  to  represent  a  single 
variable  at  just  one  atmospheric  level,  the 
40-wave  GSM  needs  3,444  terms.  This  is 
because  the  waves  are  two-dimensional- 
each  wave  requires  at  least  one  sine  and  one 
cosine  term.  For  mathematical  reasons, 


the  40-wave  model  has  41  sines  in  one 
dimension  and  42  cosines  in  the  other. 
Basically,  a  spectral  representation  finds 
the  best  mathematical  terms  using  the  sine 
and  cosine  to  describe  the  wave  motion  of 
the  observed  atmosphere. 

5.3  Resolution.  As  already  mentioned, 
the  GSM  is  a  40-wave  spectral  model  (the 
First-Guess  forecast  model  is  30-wave). 
This  means  that  the  GSM  can  recognize  40 
different  trough-ridge  combinations  in  an 
east-west  direction.  Its  application  in  the 
north-south  plane  is  much  more  difficult  to 
describe  without  the  mathematics,  but  one 
way  to  describe  the  approximate  resolution 
of  the  GSM  is  to  think  of  the  grid  on  which 
the  data  is  stored. 

The  grid  is  adapted  to  a  40’Wave  spectral 
resolution  and  is  used  for  some  of  the 
internal  calculations.  Based  on  the  grid,  its 
approximate  resolution  in  the  north-south 
direction  is  270  km;  in  the  east-west 
direction,  350  km.  The  smallest  scale 
feature  that  can  be  seen,  then,  is  about  270 
km  long  and  350  km  wide.  This  is  only 
true  for  a  surface  low  or  an  upper-level 
trough;  if  the  feature  were  a  complete 
wave  (trough-ridge  combination),  the 
smallest  feature  resolvable  would  be  twice 
that  size.  The  resolution  of  the  First-Guess 
forecast  model  is  slightly  lower. 
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TABLE  2.  Thft  Fialds  Produced  by  die  First'Guess.  Many  other  fields  are  processed  through 
algorithms  on  SDHS.  THe  list  includes  the  level^ayet^  for  which  the  data  is  available.  Ail 
the  fields  listed  are  available  for  the  00*,  06-,  and  12-hour  forecast. 


HEU>S 

directly  forecast 

Hcicht/Laycr  in  MilHbars 

Kelative  Humidity 

sfc,  1000.  ISO,  700.  SQO,  400,  300 

Heiihts 

1000,  ISO,  700,  SOO,  100,  300,  250,  200, 150, 100,  70,  50,  30,  20,  10 

Temperature 

sfc,  1000,  ISO,  700, 500, 400,  300,  250, 100,  ISO,  100,  70,  50,  30.  20,  10 

U<omponent  wind 

tfC,  1000,  ISO,  700,  SOO,  400,  300,  250,  200,  ISO,  100,  70,  SO,  30,  20,  10 

V-component  wind 

sfc,  1000,  ISO,  700,  500,  400,  300,  250,  200,  150,  100,  70,  SO,  30,  20,  10 

FIELDS 

derived 

Hei(ht/Layer  in  Millibars 

Tcmpentura  Advection 

ISO,  SDO,  300 

Spedfic  Humidity 

sfc,  1000, 150,  700,  500,  400,  300 

Tropopaiise  Heifht 

• 

Sea  Level  Pressure 

Sfc 

Dewpoint 

sfc,  1000,  ISO.  700,  SOO,  400,  300 

Vortidty  Advcoion 

1 

$00 

VortidQr 

•SO.  SOO,  300 

‘’Packed  U-V  Winds 

Sfc,  1000,  ISO,  700,  500,  400,  300,  250,  200,  150,  100,  70,  50,  30,  20,  10 

*  Tropopause  dam  varies  with  respect  to  its  exact  posirion. 

*  *  Packed  U-V  Winds  -  the  individual  U  and  V  wind  components  contain  too  many  bits  of  information  and  must  be 
compacted  for  SDHS  to  use  them. 
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•  TABLE  3.  The  Fieide  Produced  by  the  Global  Spectral  Model  (GSM).  The  first  five  are 
forecast  directly.  The  rest  are  derived  from  the  first  five.  Many  other  fields  (thickness, 
omega,  height  falls,  etc.)  aro  processed  through  algorithms  on  SDKS.  Tlie  lists  include  the 
level^ayers  for  which  ^e  data  is  available.  (*TP  =  Tropopause,  ‘‘GL  =  Gradient  Level) 
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TABLE  3,  cofitinuvd 


T/UEIUE  3,  continued. 


TABLE  3,  eonUnuftci 


TABLE  3,  coiitinu«d... 
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In  the  vertical,  the  GSM  is  a  12-layer  model. 
It  tises  a  vertical  coordinate  scheme  called 
‘'sigma  levels."  The  height  of  the  sigma 
coordinates  vary  with  the  terrain.  Uidike 
constant-pressxire  levels,  sigma  levels  never 
intersect  the  ground.  The  levels  are  closer 
together  at  about  250  mb.  The  highest  GSM 
vertical  resolution  is  between  500  mb  and 
2^  mb  (see  Figure  18)  becaxise  the  GSM  is 
an  aviation  model  and  most  aircraft  fly  at 
about  250  mb.  There  are  only  two  levels 
above  100  mb. 

The  GSM  uses  data  fields  brom  the  HIRAS 
model's  2.6*’  x  2.6®  latitude-longitude  grid  for 
its  initial  conditions.  The  HIRAS  data  is 
stored  in  the  database  on  the  GSM’s  2.5®  x 
2.5®  grid.  It  undergoes  a  spectral 
transformation  for  use  by  the  GSM  in 
making  its  forecast.  When  complete,  the 
data  *s  interpolated  back  to  the  2.5®  x  2.5® 
grid  for  storing.  Cxirrently.  most  display 
programs  require  that  the  data  be  displayed 
on  a  PST  map  projection;  this  reqtiires  the 
data  to  be  interpolated  from  the  2.5®  X  2.5® 
degree  grid  to  file  whole-mesh  subgrid  and 
results  in  a  further  decrease  in  data 
resolution.  Products  produced  at  AFGWC 
on  SDHS  (or  for  transmission  over  AWDS  or 
facsimile  circuits)  are  based  on  the  lower 
resolution  whole-mesh  subgrid. 

5.4  Normal  Mode  Initialization.  Before 
the  GSM  can  use  the  HIRAS  analysis  data 
for  its  initial  conditions,  it  must  undergo  a 
process  called  "Normal  Mode  Initialization," 
or  NMI.  The  main  purpose  of  the  NMI  is  to 
prevent  the  model  from  becoming 
mathematically  unstable  due  to  inherent 
instabilities  in  the  input  data  that  are  a 
result  of  minor  perturbations  the  model 
can't  resolve.  These  perturbations  may  be 
real  (small-scale  phenomena),  fictitious,  or 
due  to  small  observation  errors.  Almost  all 
models  use  some  form  of  NMI. 
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Figure  18.  The  Sigma  Levels  of  the  Globai 
Spectral  Model.  This  is  the  vertical 
resolution  of  the  model.  Note  the  higher 
concentration  of  layers  between  600  and  200 
mb.  The  forecast  of  some  elements  is 
conducted  directly  on  the  layer  boundaries 
(indicated  by  the  thick  horizontal  lines). 
The  rest  of  the  forecast  is  conducted  in  the 
center  of  mass  between  the  layer  boundaries 
(from  Stobie  1986). 
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S.S  Smaii-Scaie  Dlffuslomi.  This 
mechanism  is  used  duriu^  the  running  of 
the  model;  it  also  is  intended  to  keep  it 
firom  becoming  mathematically  unstable. 
Diffusion  dampens  the  higher  frequency 
horizontal  waves  (e.g.,  gravity  waves)  ihat 
the  model  cant  resolve.  A  simplistic  way  to 
put  it  is  that  it  smooths  the  data  fields  very 
slightly.  In  the  GSM,  the  need  for  small- 
scale  diffiudon  is  not  great,  since  the  GSM 
does  not  t^cle  back  on  itself-that  is,  the 
GSM  does  not  use  its  12-hour  forecast  to 
initialize  the  next  forecast  run.  The 
First-Guess  model,  on  the  other  hand,  does 
cycle  back  on  itself  and  therefore  requires 
much  more  diffusion. 


5.6  Topography.  The  topography 
incorporated  in  the  GSM  is  much  smoother 
than  actual  terrain  (see  Figure  19).  There 
are  no  jagged  peaks  nor  deep  valleys,  and 
the  Earth  appears  to  be  composed  of  rolling 
hills.  Again,  the  main  reason  for  using  the 
smoother  topography  is  to  prevent  the  model 
from  becoming  numerically  unstable.  Real 
terrain  data  would  create  features  within 
the  forecast  that  the  model  could  not 
resolve.  Another  reason  for  the  use  of  the 
smoothed  terrain  is  to  optimize  the 
forecasting  capability  of  the  computer,  given 
the  limited  memory  and  CPU  available. 


Figuro  19.  The  Topography  of  North  America  as  Seen  by  the  GSM.  The  heights  are  in 
meters  (from  Stobie,  1986), 


5.7  Surface  Friction.  Becatise  of  surface 
frictiion,  winds  in  the  GSM  are  gfenerally 
nou'geostropbic.  Over  major  mountain 
chains  (such  as  the  Rockies  and  Himalayas) 
suarface  friction  is  considered  to  be  five  times 
stronger  tiian  over  oceans. 

5.8  Sensiblo  Heat  Exchange.  In  the 

GSM,  the  atmosphere  receives  sensible  heat 
only  from  the  oceans.  No  heat  is  exchanged 
over  land  surfaces.  Sensible  heat  released 
into  the  atmosphere  is  dependent  on 
temperature  differences  between  the  water 
and  the  air,  and  on  the  surface  wind  speed. 
If  the  water  is  warmer  than  the  air,  heat 
energy  is  transferred  from  the  water  to  the 
air,  and  vice-versa.  When  strong  winds  are 
forecast,  the  sensible  heat  exchange  is  two 
to  three  times  higher  than  during  calm 
conditions.  The  heat  received  or  lost  by  the 
atmosphere  is  based  on  monthly 
climatological  sea-surface  temperatures. 
When  an  £1  Nino  (or  another  dramatic 
anomaly  in  sea-surface  temperattue)  occurs, 
the  (^M  will  not  accurately  represent  the 
sensible  heat  exchanges;  the  result  is  an 
occasional  forecast  error  in  cyclone 
strengthening  or  weakenmg.  The  GSM 
would  only  recognize  the  heat  already 
exchanged  into  the  atmosphere  as  evidenced 
by  the  analysis. 

5.9  Evaporation.  Only  moisture  from  the 
oceans  can  be  evaporated;  the  GSM  does 
not  take  soil  moisture  evaporation  into 
account.  Evaporation  is  based  on  monthly 


climatologicail  sea-surface  temperatures  and 
forecast  surface  wind  speeds.  Stronger 
winds  result  in  increased  evaporation. 
Because  of  these  limitations,  GSM  moisture 
forecasts  are  typically  too  dry,  especially 
over  land. 

5.10  Precipitation  Effects.  The  GSM 

forecasts  two  types  of  precipitation:  synoptic 
(large-scale)  and  cumulus  (small  scale). 
Evaporation  and  the  release  of  latent  heat 
residtmg  from  precipitation  are  taken  into 
account  during  the  model  forecast  process. 
The  GSM  is  not  considered  a  superior 
precipitation  model  and  these  processes 
don't  help  the  moisture  forecasts.  The  real 
purpose  of  GSM  precipitation  forecasting  is 
to  forecast  the  effects  of  precipitation  on 
other  forecast  elements;  that  is,  to  quantify 
latent  heat  release  and  evaporation  and 
determine  its  effects  on  temperature  and 
wind. 

5.11  Radiation.  This  energy  source 
currently  does  not  exist  in  the  GSM;  there 
is  no  radiation  cooling  at  night  nor 
imolation  during  the  day.  Although  the  net 
result  is  a  warm  bias  in  the  model,  it  is  not 
considered  overly  signiAcant  in  the  short 
range  (3  days  or  less).  The  First-Guess 
forecast  model  is  much  more  vulnerable  to 
the  lack  of  radiation  processes  because  it 
cycles  back  on  itself.  The  problem  is  kept 
under  control  through  Optimum 
Interpolation  by  blending  the  First-Guess 
forecasts  with  observations  every  6  hours. 


5.12  GSM  Performance.  Not  much  has 
been  written  on  the  GSM's  capabilities  and 
weaknesses.  Most  of  the  data  for  the 
following  notes  on  GSM  performance  was 
collected  firom  statistical  processes,  but  some 
of  the  information  is  based  on  discussions 
with  AFGWC^YS  meteorologists. 

•  Long  waves  are  handled  well,  considering 
that  the  model  was  designed  for  flightdevel 
forecasting. 

•  Forecast  wind  speeds  are  typically  too 
slow  below  the  100-mb  level. 


•  Moisture  forecasts  are  typically  too  dry. 

•  The  GSM  is  generally  too  warm  in  the 
troposphere. 

•  Since  the  GSM  forecasts  geopotential 
heights  balanced  with  temperature,  the 
geopotential  heights  are  generally  too  high. 
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Chapter  6 

THE  REAL-TIME  NEPHANALYSIS  MODEL  (RTNEPH) 


6.1  introductioil.  The  main  purpose  of 
the  RTNnEiFH  model  is  to  initialize  the  cloud 
forecast  models.  E^TNEPH  uses  an 
eighth-mesh  grid.  Its  horizontal  resolution 
is  about  48  km  at  60°  N.  RTNEPH  has  only 
four  levels  in  the  vertical.  These  levels  are 
not  fixed,  but  are  sorted  by  cloud  base;  the 
highest  base  is  at  level  one.  The  four  levels 
are  not  necessarily  continuous  horizontally, 
and  may  be  at  different  levels  for  each 
gridpoint.  If  there  are  more  than  four  cloud 
layers,  extra  layers  are  merged  with  the 
next  closest.  Up  to  eight  layers  can  be 
merged  to  fit  into  the  four  available  levels 
per  grid  point.  If  there  are  more  than  eight 
layers,  the  lowest  wiU  not  be  accoimted  for. 
Levels  can  range  from,  the  surface  to  21,900 
meters  (71,932  feet)  MSL.  The  vertical 
precision  can  be  as  fine  as  30  meters 
between  levels  below  6,000  meters  and  300 
meters  between  each  level  above  6,000 
meters. 

6.2  RTNEPH  Production  Cyclea.  There 
are  three  RTNEPH  production  cycles: 
sprint,  non-sprint,  and  update.  The  sprint 
cycle  is  only  used  when  support 
requirements  demand  it.  The  other  two  are 
run  every  day.  Unlike  HIRAS,  the 
RTNEPH  data  base  is  updated  continuously 
using  these  three  cycles. 

0.2.1  The  Sprint  Cycle.  The  purpose  of  the 
sprint  cycle  is  to  take  a  quarter  orbit  of 
satellite  data  as  soon  as  it  is  received  and 
incorporate  it  into  the  RTNEPH  analysis  for 
immediate  tise  in  the  short-range  cloud 
forecast  model.  The  sprint  cycle  is  quality- 
controlled  manually  as  needed.  Changes  are 
made  with  bogus  data.  Bogus  data  is 
information  manually  given  to  the  analysis 


model  to  counter  the  negative  affects  of  bad 
conventional  or  satellite  data.  The  sprint 
cycle  only  updates  that  portion  of  the 
RTNEPH  analysis  for  which  satellite  data  is 
available,  but  both  satellite  and 
conventional  data  are  used  to  update  the 
RTNEPH  for  these  points. 

6.2J!  The  Non-Sprint  Cycle.  The  non-sprint 
cycle  is  the  same  as  the  sprint  cycle  mccept 
for  the  time  factor.  Immediate  cloud 
forecast  model  input  and  manual  quality 
control  are  not  needed.  The  non-sprint  cycle 
is  still  triggered  by  the  receipt  of  a  quarter- 
orbit  of  satellite  (kita.  Like  the  sprint  cycle, 
it  only  updates  that  portion  of  the  RTNEPH 
analysis  for  which  satellite  date  is  available. 
The  non-sprint  cycle  also  uses  satellite  and 
conventional  data  to  update  the  RTNEPH. 

641.3  The  Update  Cycle.  The  update  cycle  is 
run  every  3  hours  to  incorporate  as  much 
data  into  the  RTNEPH  database  as  possible. 
Unlike  the  other  two  cycles,  the  update  cycle 
updates  the  entire  Northern  and  SouUiem 
Hemispheres  during  its  trun.  The  Northern 
Hemisphere  is  updated  1  1/2  hours  after 
data  time;  the  Southern  Hemisphere,  2  1/2 
hours  after  data  time.  Conventional  data  is 
the  primary  input  because  the  satellite  data 
(received  in  quarterly  segments)  has,  for  the 
most  part,  already  been  used  to  update  the 
analysis.  All  conventional  data  received 
from  both  hemispheres  is  incorporated;  in 
the  other  two  cycles,  conventional  data  is 
only  used  for  the  quarter  orbit  being 
processed. 

6.3  Analysis  Modes.  Theire  are  three 
types  of  RTNEPH  analyses:  conventional, 
satellite,  and  persistence. 


9  J.1  ConvMitkMtid  AiuUyal*.  Since  surface 
olMMarvations  are  updated  in  the  database 
every  hour,  the  conventional  cloud  analysis 
is  a^  perfonued  once  an  hour.  This  allows 
for  tike  miMt  recent  data  to  be  available 
whenever  tiie  Merg0  Processor  (to  be 
explained  later)  needs  it.  Conventional  data 
is  only  used  at  those  gridpoints  for  which  it 
is  avakilable.  OX  is  not  performed. 

Conventional  cloud  data  from  upper-air  and 
surface  synoptic,  MjEHi’All,  and  AIRWAYS 
reports  is  used  (including  the  data  unique  to 
each  type  of  observation).  If  more  than  one 
source  of  conventional  data  is  available, 
conflicts  are  resolved  by  taking  into  account 
the  data's  timeliness  and  reliability. 
Surface-baaed  observations  are  considered 
more  reliable  than  upper-air  soundings  for 
the  same  point.  All  hourlies  and  special 
observations  3  hoturs  prior  to  analysis  are 
oknsidered. 

When  doud  amokints  are  reported  as 
scattered,  broken,  or  overcast  in  surface 
reports,  the  data  is  stored  directly  and 
assigned  25,  76,  and  1(X)%  coverage  values, 
respectively.  Cloud  bases  reported  in 
Airways  code  are  converted  to  meters. 
Because  doud  information  in  synoptic 
observations  is  limited,  cloud  height 
eateries  are  used  (i.e.,  surface  to  6,5(X) 
feet,  6,500  to  22,000  feet,  and  22,000  feet 
and  above).  If  a  doud  type  is  not  observed, 
zero  percent  probability  is  assigned  to  the 
corresponding  height  category.  When  a 
valid  doud  type  is  given,  1(X)%  probability  is 
assigned.  For  missing  doud  types,  50% 
probability  is  assigned. 

Present  weather  is  also  vised  to  determine 
doud  base  heights  and  types.  Only  one  type 
of  predpitation  (the  most  significant)  can  be 
included  in  the  RTNEPH  database  for  each 
gridpoint.  It  is  maintained  in  the  database 
as  a  separate  element. 


Model  terrain  is  used  to  change  AGL  to 
MSL;  this  only  affects  data  for  individual 
gridpoints. 

6.3J1  Sateilita  Anaiyais.  Satellite  information 
is  only  used  when  polar  orbiting  data 
(infrared  an<Vbr  visible)  is  available.  The 
threshold  method  for  determining  doud 
cover  from  satellite  data  is  the  basic  method 
used.  Eackgrovmd  brightness  fields  are  used 
as  a  starting  point  for  a  histogram.  The 
ability  to  tell  the  differ^ce  between  doud 
cover  and  background  is  extremely 
dependent  on  these  fields,  which  change 
with  sun  angle  and  reflectivity  (snow/ice) 
and  must  be  updated.  Sensed  brightness  of 
a  pixel  determined  to  be  dear  is  compared 
with  the  background  brightness,  which  is 
adjusted  so  that  it  approaches  the  value  of 
the  new  sensed  value  without  straying  too 
far  from  the  current  background  brightness 
values. 

In  general,  darker  Cess  bright)  visual  data 
is  associated  with  cloud-firee  areas. 
Conversely,  brighter  data  is  associated  with 
doudy  areas,  but  brighter  values  may  also 
be  associated  with  terrain  albedo  affects,  or 
snow  and  ice.  Infirared  data  represents 
temperature  in  terms  of  brightness.  In 
RTN£PH,  brighter  infrared  measuramenta 
indicate  lower  temperatures.  Although  both 
visual  and  infrared  provide  brightness 
values  indicating  amounts  of  reflected  or 
radiated  energy,  the  determination  of  doud 
cover  requires  a  comparison  of  the  satellite 
data  to  the  expected  background  brightness. 
For  visual  data,  the  background  is  a 
variation  of  surface  albedo.  Therefore,  if  the 
sensed  data  is  brighter  than  the  background, 
clouds  are  present.  For  infrared  data, 
backgrovmd  brightness  is  based  on  surface 
temperature;  if  the  sensed  data  is  colder 
than  the  adjusted  backgrovmd  surface 
temperatures,  clouds  are  present. 
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When  visual  satellite  data  is  received  by 
HTNEFH,  the  pixel  data  would  appear  as 
shown  in  Step  2  of  Figure  20.  Step  3  shows 
an  array  of  the  visual  bi^htness  values  for 
the  pixel  gray  shades  from  the  Satellite 
Global  Database  (SGDB).  Step  4  shows  the 
resultant  histogram,  which  gives  frequency 
of  occurrence  for  the  possible  gray  shades. 


After  the  histogram  is  formed,  the  expected 
background  brightness,  (the  value  of  5 
repi*esents  a  water  background  in  this 
example)  is  compared  to  the  histogram.  All 
pixels  with  a  value  highei^brighter  than  the 
expected  background  are  considered  "cloudy" 
pixels. 
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Figure  20.  A  Pictorial  Example  of  How  Cloudy  Pixels  are  Determined  by  the  Satellite 
Processor  (from  Klees  and  Cox,  1988). 
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BJ3A  PtrstotniM  Aiurfysto.  The  STNEPH 
p4»rsistence  analysis  is  simply  one  that  is 
equal  to  the  previous  analysis.  The 
persistence  analysis  is  used  when  no  new 
data  is  available. 

6.4  Th*  M«rge  Proceasor.  Ihe  three 
analj^ses  described  above  are  computed  in 
three  separate  processors.  When 
compl^«d,  they  are  blended  together  in  the 
merge  proc^sor.  Since  conventional  data  is 
input  only  at  the  gridpoints  for  which  they 
are  available,  and  since  satellite  data  is 
input  only  where  it  is  available,  there  will 
be  some  areas  with  no  data;  these  are 
where  the  persistence  data  will  be  input. 
When  the  merge  processor  is  finished,  it 
provides  an  updated  global  nephanalysis 
database.  The  merge  processor  also 
reconciles  conflicts  between  conventional 
and  satellite  data  at  the  same  grid  point. 

6.4.1  Bogus  Data.  If  bogoa  data  is  available 
for  a  particular  gridpoint,  it  will  replace  all 
other  data  only  if  timeline  requirements  are 
met. 

6.4.2  Dida  Spreading.  The  merge  processor 
also  considers  conventional  data  that  can  be 
spread  to  other  gridpoints.  Data  is  spread 
fix>m  gridpoints  that  are  flagged  as  "best 
report"  gridpoints  using  the  following 
guidelines: 

•  Data  will  be  spread  only  to  those 
gridpoints  within  an  allowable  spread 
radius. 

•  Data  is  spread  firom  the  nearest,  "best 
report"  gridpoint. 

•  If  the  distances  between  two  "best  report" 
gridpoints  are  equal,  the  newest  report  is 


used.  If  they  have  the  same  time  value,  the 
least  cloudy  report  is  used.  If  they  axe 
equally  cloudy,  the  one  with  the  greatest 
visibility  is  used.  If  no  tie  breaker  is  found, 
then  the  first  "best  report"  considered  is 
used. 

6.5  RTNEPH  Deficiencies. 

•  When  only  satellite  data  is  available, 
RTNERI  under-interprets  low  clouds 
because  the  infirared  threshold  becomes 
laiger  than  the  temperature  diff^ence 
between  the  clouds  and  the  ground  or 
between  the  ground  and  the  modeled  ground 
temperature. 

•  Over-  or  under-interpretation  may  occur 
along  coastlines  due  to  problems  in  choosing 
a  representative  background  brightness  or 
temperature  to  represent  the  gridpoint 
during  the  satellite  analysis.  For  example, 
sea  and  land  surface  temperatures  are 
almost  always  different.  The  satellite 
processor  building  the  histogram  may  use 
higher  daytime  land  surface  temperatures 
as  a  background  over  the  water  and 
therefore  see  clouds  where  there  aren't  any. 

•  Snow  and  ice  fields  are  incorporated  into 
RTNEPH  by  the  SNODEP  model.  If  the 
snow  and  ice  data  is  incorrect,  the  RTNEPH 
might  analyze  stratus  instead  of  snow. 

•  High  plateaus  (such  as  the  Tibetan 
Hateau)  are  generally  over-interpreted 
because  of  their  lower  surface  temperatures. 

•  Small-scale  clouds  are  usually 
under-interpreted  due  to  the  satellite  pixel 
resolution. 
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Chapter  7 

THE  SLAYER  MODEL 


7.1  Model  Description.  The  SLAYER 
Model  produces  forecasts  of  clouds, 
temperatures,  and  other  elements  at  3>hour 
time  steps  out  to  48  hours  in  the  Northern 
Hemisphere  and  out  to  24  hours  in  the 
Southern  Hemisphere.  5LAYER  produces 
f<arecasts  on  an  octagon  half*mesh  grid.  It 
does  not  completely  cover  the  tropics. 
SLAYER  forecasts  for  the  following 
atmospheric  layers: 

•  The  gradient'level  layer  (600  to  987 
meters  AGL) 

•  The  850*mb  layer  (^7  to  2,235  meters 
MSL) 

•  The  700-mb  layer  (2,235  to  4,238  meters 
MSL) 

•  The  5(X)-mb  layer  (4,239  to  7,368  meters) 

•  The  300>mb  layer  (7,368  meters  MSL 
and  above). 

Before  running  a  forecast,  5LAYER 
compacts  (or  smooths)  the  eighth-mesh 
RTNEFH  analysis  into  a  half-mesh,  lOO-NM 
resolution  grid  analysis.  RTNEPH,  GSM, 
and  HIRAS  are  the  primary  data  sources  for 
SLAYER,  which  uses  dew-point  depression, 
static  stability,  and  cloud  type  to  produce 
forecasts  of  clouds,  precipitation  amount, 
precipitation  type,  and  icing. 

7.2  Resolution.  As  its  name  implies, 
5LAYER  forecasts  for  five  layers  in  the 
vertical  on  a  half-mesh  grid  with  a 
horizontal  resolution  of  190.5  km  at  GO**  N. 
Although  the  resolution  of  SLAYER 
forecasts  is  about  equal  to  that  of  many 
operationsd  dynamic  models,  the  forecasts 
are  smoothed  considerably  compared  to  raw 


satellite  data,  which  resolves  cloud  elements 
to  3  km.  The  lower  SLAYER  resolution 
results  in  a  final  forecast  that  represents 
average  cloud  conditions  within  a  la3rge 
volume.  In  later  discussions,  this  volume  is 
referred  to  as  a  "parcel."  Vertical 
atmospheric  layers  (850  mb,  700  mb,  500 
mb,  and  3(X)  mb)  are  fixed,  while  the 
gradient  level  follows  terrain.  Fixed  levels 
may  intersect  the  ground  at  higher 
elevations  where  the  gradient  layer  is  above 
850  or  700  mb. 

7.3  The  Quasi-Lagrangian  Advection 
Scheme.  The  model  makes  its  forecast  by 
advecUng  parcels  (or  volmnes)  of  the 
atmosphere  using  a  quasi-Lagrangian 
advection  scheme-essentially  a  trajectoxy 
forecast.  A  tiue  Lagrangian  scheme  takes  a 
parcel  firom  the  initial  point  to  its 
destination  at  the  end  of  the  entire  forecast 
period.  The  SLAYER  takes  a  parcel  from  its 
initial  point  to  the  first  3-hour  time  step, 
then  to  the  next  3-hour  time  step,  and  so  on. 
The  KTNEIPH  provides  the  initial  cloud 
information  for  each  gridpoint.  The  GSM 
provides  the  winds  for  the  initial  time  and 
for  each  S-hour  time  step,  and  HIRAS 
provides  the  initial  temperature  and 
moisture  conditions.  5LAYER  modifies  air 
parcels  as  they  move  through  the  trajectory 
path  to  the  3-hour  point.  The  3-hour 
forecast  then  becomes  the  initial  state  for 
the  next  3-hour  forecast.  This  process  is 
repeated  until  the  desired  forecast  length 
(48  houi^)  is  reached. 

7.4  Initializing  the  SLAYER.  The 

RTNEPH  analysis  provides  cloud  top,  base, 
and  amoimt  information  for  each  layer  at 
each  gridpoint.  Total  cloud  coverage  (in 
percent)  is  also  provided.  In  order  to 
initialize  5LAYER  using  RTNEPH  data,  the 


four  xum>fixed  levels  of  the  BTNEFH 
anal3mis  data  must  be  matched  to  the  five 
layrans  tibee  forecast  model.  The  data  must 
also  be  cmupacted  (firom  the  BTNBPH 
aaa^ysis  eighth-mesh  grid  to  the  51AYER 
half-iaeah  gricQ  to  reduce  the  volume  of  data 
while  retainmg  most  of  the  larger-scale 
infiarmation. 

Ihe  compacfion  takes  data  from  the  25 
e^th-mesh  gri^oints  axid  applies  it  to  the 
oae  half-mesh  gridpoint  they  surroimd.  It 
uses  a  grid>weighfi2ig  system  in  which  all 
points  that  immediately  surround  the 
gridpoint  on  the  SLAYER  have  equal 
weight;  eighth-mesh  gridpoints  farther 
away  from  the  6LAYER  half-mesh  gridpoint 
affect  it  to  a  lesser  degree.  This  proems 
degrades  the  analysis  because  data 
averaging  results  in  a  loss  of  distinction 
between  moist  and  dry  air.  To  spread  the 
clouds  vertically,  the  top  and  base  of  each 
BTNEFH  cloud  layer  is  checked  to 
determine  which  5LAYER  fixed  layers  they 
overlap.  BTNEFH  cloud  layer  amounts  are 
inserted  into  5LAYER  fixed  layer  gridpoints 
unless  the  receiving  gridpoints  already  have 
a  larger  cloud  amount.  For  5LAYER  to  use 
the  eighth-meah  BTNEFH  analysis,  it  would 
need  to  initialize  2,621,440  grid  points 
(262,144  grid  point^hemisphere  times  two 
hemispheres  times  five  layers).  Computer 
hardware  and  AFGWC  time  constraints 
make  it  impossible  for  a  prediction  model  to 
use  eighth-mesh  gridpoints  at  this  time. 

7.5  Old  D«Ul  Limitations.  In  areas 
where  the  BTNEFH  is  too  old  (by  4  ho\irs  or 
more),  5LAYER  uses  its  previous  forecast  as 
an  initial  condition,  or  analysis.  This 
commonly  occurs  in  extratropical  areas 
where  polar  orbiting  satellite  data  can  be  5 
1/2  hours  old,  or  more.  Fast-moving  cloud 
systems  result  in  a  loss  of  distinction 
between  cloudy  and  dry  air  due  to  the 
persistence  forecast  initialization  process  in 
areas  of  old  data. 


7.6  Initializina  the  RTNEPH.  Total  cloud 
amount  for  a  particular  vertical  grid  stack  is 
computed  h*om  the  statistical  union  of  layer 
cloud  amotmts,  maximum  layer  cloud 
amount,  and  average  layer  separation.  To 
help  understand  this,  let's  assume  that 
there  are  two  layers,  both  with  50%  cloud 
coveraj^.  If  they  are  separated  by  11,000 
meters  (assumed  to  be  the  depth  of  the 
troposphere)  the  layers  are  considered 
uncomlated-Uxe  total  doud  amount  for  the 
stack  is  97%.  If  the  two  layers  are 
immediately  adjacent  to  one  another,  they 
are  considered  correlated-the  total  doud 
amoiint  is  equal  to  the  maximum  doud 
amount  of  the  individual  layers,  or  50%. 

If  11,000-meter  separation  of  layers  provides 
a  97%  total  cloud  amount,  and  if 
immediately  adjacent  layers  result  in  doud 
amounts  equal  to  the  maximtun  doud 
amount  between  the  two  layers  in  percent, 
layeia  with  separations  between  these  two 
figures  are  computed  proportionally-that  is, 
50%  maximum  cloud  with  two  layers; 
6,500'meter  separation  equals  about  74% 
total  coverage. 

Where  more  than  two  layejrs  are  found,  the 
average  of  all  possible  pair-vtise  separations 
of  the  doud  layers  is  divided  by  11,000,  then 
subtracted  from  1  to  derive  the  percent  of 
total  cloud;  1  is  equal  to  100%.  Hub  new 
cloud  amount  is  used  to  adjiut  the  5LAYER 
amounts  to  match  reality  (the  BTNEFH).  If 
the  total  doud  amount  is  not  close  enough  to 
the  analyzed  total  as  compared  to  BTNEFH, 
they  are  considered  to  be  out  of  tolerance 
and  corrected  with  a  computation  that 
includes  the  BTNEFH  total  cloud  amo\mt. 

7.7  Condensation  Pressure  Spread 
(CPS)  is  the  moisture  variable  used  by  the 
SLAYER  CPS  is  the  difference  in  pressure 
(millibars)  between  a  parcel  and  the 
pressure  to  which  it  needs  to  be  lifted  for 
condensation  to  occur.  SLAYER  assvunes 
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iiiat  if  a  parcel  is  lifted  to  a  low  enotigh 
pressure,  ctmi^nsation  will  occur.  Cl^ 
oMTelates  directly  to  dew-point  depression. 
Altkough  dew  points  are  not  advected, 
I»a!vel8  are  moved  in  a  trajectory,  dew 
points  are  depmident  on  the  pressure  of  ^e 
parcel  at  tbke  end  of  the  forecasted  move. 
Since  tbe  BTNEPH  cant  determine  the 
moisture  content  of  cloud-ftee  areas,  either 
the  default  CPS  or  Uie  CPS  derived  from  the 
GSM  (whichever  is  drior)  is  used  to  fill  tiie 
gaps  around  cloudy  areas.  CPS  is  also  the 
link  between  the  cloud  amount  data  received 
from  PTNEPH  and  the  GSM  dew-point 
data,  making  it  a  "common  to  both" 
moisture  variable. 

For  example,  a  CPS  of  5  would  equal  97% 
cloud  coverage  and  a  dew-point  depression 
of  zero.  An  advantage  to  having  a  CPS  is 
that  it  correlates  changes  in  cloud  amount 
with  changes  in  vertical  motion.  If  an 
unsaturated  parcel  is  lifted,  adiabatic 
cooling  takes  place,  the  dew-point  depression 
narrows,  and  clouds  form.  For  downward 
motion,  the  dew-point  depression  widens 
and  clouds  dissipate. 

Tlie  cloud  amount  can  be  determined 
directly  by  the  net  vertical  displacement  of 
the  parcel  in  terms  of  pressure.  If 
supersaturation  occurs  due  to  the 
modification  of  a  parcel  along  its  trajectory, 
the  excess  moisture  is  condensed  and  added 
to  the  quantitative  precipitation. 

7.8  Cloud  Type  Determination. 

Determination  of  cloud  type  in  the  5LAYER 
is  a  relatively  simple  process.  Cumuliform 
or  stratiform  type  clouds  are  determined  by 
the  stability  of  the  850-,  700-,  and  500-mb 
layers.  F^or  example,  if  the  wet-bulb 
potential  temperature  of  the  700-mb  layer  is 
less  than  that  of  the  850-mb  layer,  positive 
(upward)  buoyancy  in  the  700*mb  layer  is 
assumed. 


When  there  is  positive  buoyancy  in  a  layer 
containing  clouds,  the  cloud  type  is 
cumuliform;  if  there  is  negative  buoyancy, 
the  cloud  type  is  stratiform.  Lf,  in  a  vertical 
column,  there  are  clouds  in  a  xwgatively 
buoyant  700-mb  layer,  the  cloud  type  is 
altostratus;  if  the  700-mb  layer  is  positively 
buoyant,  the  cloud  type  is  altocumulus. 

If  the  850-,  7(X}-,  and  500-mb  lay^  are 
positively  buoyant  end  all  three  contain 
clouds,  the  cloud  type  is  cumulonimbus.  If 
the  column  is  negatively  buoyant,  cloud  type 
is  nimbostratus.  Cimis-level  douds  are 
considered  simply  "cirrus,"  not  cirrostratus 
or  cirroctuxiulus. 

7.9  Precipitation  Type  Determination. 

Cloud  type,  as  determined  lAW  7.8,  is 
primarily  used  in  determining  precipitation 
type.  Other  variables  used  are  amount  of 
cloud  cover  and  static  stability. 

Showery  precipitation  requires  the  presence 
of  towering  cumulus  or  cumulonimbus,  and 
at  least  25%  doud  cover.  To  forecast 
snowshowers,  the  gradient  level  layer 
temperature  must  be  0°  C.  The  intensity  of 
showery  precipitation  is  determined  from 
the  static  stability,  which  is  determined  by 
computing  the  Showalter  stability  index  for 
gridpoints  with  terrain  heights  below  850 
mb.  For  terrain  above  850  mb,  stability  is 
detennined  in  the  same  way  as  for  doud 
type.  In  areas  of  cumulonimbus,  the 
l^otal-Totals  stability  index  is  used. 

Diurnal  heating  and  orographic  effects  are 
not  considered  in  determining  precipitation 
type.  Continuous  precipitation  is  forecast  if 
cloud  cover  is  at  least  75%  and  there  are 
cloud  types  other  than  cumulonimbus  or 
towering  ciimulus  in  the  860-  and  700-mb 
layejrs. 
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'Hja  saawAndxii  di8tinctio&  i«  made  in  the 
same  way  as  for  showery  precipitation. 
IntKudty  dbterm£ctati<m  is  a  function  of  total 
moisture  depth  and  the  buoya33£y  of  the  700- 
mh  layex^  the  greater  the  total  moisture 
d^th  and  700-mb  buoyancy,  the  greater  the 
precipitatiim  intensity. 

7.10  TempemUiiw  InitiiUizatfon.  The 
initial  temperature  field  needed  for  the 
SLAYER  forecast  is  input  from  the  KIBdJS 
anal$^  and  converted  frtun  the  HIRAS  2.5** 
X  2.5°  degree  grid  to  a  whole-mesh  grid, 
which  is  then  converted  to  the  half-mesh 
grid  used  by  the  SLAYER.  Hm 
temperatures  are  interpolated  direcUy  to 
ea(h  the  four  pressure  levels  that  define 
the  midpoint  of  Ihe  fixed  model  layers. 

The  terrain-following  gradient  level  is 
initialized  using:  a  height  interpolation 
between  fbmd  preesure  levels  in  HIRAS  to 
ensure  that  gradient-layer  temperatures  are 
relatively  free  from  amall-scale  local  effects. 

A  trajectory  method  of  temperature 
advection  is  used  to  make  temperature 
forecasts.  As  tngectories  displace  parcels, 
they  are  modified  by  several  atmospheric 
processes-the  dry  adiabatic  process  is  tiie 
most  important  of  these,  but  others  include 
latent  heat  release,  diurnal  heating,  and 
turbulent  mixing.  Some  of  these  processes 
are  not  handled  by  the  GSM.  The 
temperature  lapse  rate  at  each  gridpoint  is 
checked  at  the  end  of  each  time  step  to 
make  sure  a  superadiabatic  lapse  rate  was 
not  created. 

7.11  Fomeast  Etemants,  Temperature 
and  CPS  are  the  only  two  meteorological 
elements  advected  using  the  trajectory 
scheme.  Other  fields  derived  from  these  two 
elements  include  layered  and  total  cloud 
cover,  dew-point  depression,  stability, 
quantitative  precipitation,  cloud  type,  and 
icing. 


7.12  Wind  initiiiailEation.  Wind  forecasts 
from  the  GSM  are  tued  to  4»nstruct  three- 
dimensional  trajectories  to  determine  the 
path  of  each  air  parcel.  After  the  wind 
forecasts  are  averaged  in  time,  the 
trajectories  are  computed  backward  firom 
selected  gridpoints  to  their  origins.  Since 
the  trajectori^  define  the  origin  of  an  air 
parcel,  they  are  <K>mputed  by  starting  at  the 
final  forecast  point.  Ibis  allows  for 
modification  of  the  parcels  through  time  and 
space. 

7.13  Terrain  effects  are  determined  using 
the  half-mesh  terrain  database.  Since  the 
gradient  level  is  terrain  following,  the  wind 
components  received  firom  the  GSM  may  be 
mof^ed  to  account  for  steep  terrain,  for  two 
reasons: 

•  First,  to  prevent  an  upstream  parcel 
trajectory  firom  backing  into  the  terrain.  At 
the  conclusion  of  each  time  step,  the 
pressure  of  the  trajectory  origin  is  compared 
to  the  interpolated  tenrain  pressure.  If  the 
origin  intersects  (or  is  below  the  terrain 
level),  the  vertical  component  of  the 
trajectory  is  adjusted  so  that  it  is  exactly  20 
mb  above  the  terrain. 

•  Second,  when  a  large  change  in  pressure 
occurs  in  a  short  time  (as  when  air  parcels 
move  over  mountains),  the  model  has  a 
tendency  to  create  heavy  precipitation  on 
the  windward  side  and  counter  with  extreme 
dryness  on  the  leeward  side.  The  vertical 
displacement  of  each  gridpoint  is  checked  to 
see  if  it  exceeds  the  50-mV3-hour 
displacement  criterion;  if  i^  does,  the 
horizontal  component  of  the  tra  .;tory  is  cut 
by  one  half  and  the  vertical  displacement  is 
recomputed.  If  the  criterion  is  still 
exceeded,  the  horizontal  component  is  again 
cut  by  one  half  and  the  vertical 
displacement  is  computed  again.  Ibe 
process  is  repeated  until  the  displacement 
value  is  less  than  50  mV3  hours. 


45 


7.14  The  SLAYER  Grid.  The  grid  used 
for  dre  SLAYER  is  the  octegozual  half-mesh 
lhat  was  shown  in  Figure  5.  This  poses  a 
problem  when  trajectory  origins  are  beyond 

boundary  of  the  grid.  The  SLAYER 
sol'mi  this  problem  by  truncating  all 
trajectories  for  the  sides  of  the  grid  that 
parallel  the  x-y  axis  {x  is  170“  E,  10“  W)  to 
dxe  edge  of  the  octagon  grid.  More 
complicated  mathematics  axe  applied  for  the 
diagonal  sides.  Simplistically  put,  the 
diagonal  sides  are  rotated  45“,  the 
truncation  is  completed,  and  the  sides  are 
rotated  back.  All  trajectories  within  two 
whole-mesh  grid  points  (400  NM)  of  the 
edge  ixiside  the  grid  botmdary  are  checked  to 
see  if  this  truncation  is  needed.  No  parcel 
irom  outside  the  boundary  of  the  grid  is 
allowed  to  move  across  the  boundary  into 
the  grid,  and  flow  is  forced  to  be  parallel  to 
the  boundary  of  the  grid. 

7.15  Trajectory  Origins.  During  the 
forecast  process,  trajectory  origins  rarely  fall 
directly  on  a  gridpoint  at  one  of  the  fixed 
layers  of  the  model.  Because  of  this,  two 
forms  of  data  interpolation  have  been 
developed:  two-point  and  eight-point.  Using 
both  methods  gives  a  precise  description  of 
areas  with  large  amounts  of  cloud  coverage 
and  decreases  computer  processing  time. 

•  Two-point  interpolation  is  xwed  when  the 
parcel  originated  hxtm  an  area  where  the 
analysis  has  50%  cloud  cover  or  less.  In  the 
two-point  interpolation,  a  parcel's  point  of 
origin  is  rounded  to  the  nearest  vertical 
gridpoint  axis.  Data  is  then  interpolated  to 
the  point  from  the  nearest  two  hxed  layers. 
Because  of  rounding,  two-point  interpolation 
tends  to  be  less  accurate  than  eight-point 
interpolation. 

•  Eight-point  interpolation  is  used  when  the 
analysis  shows  more  than  50%  cloud 
coverage.  It  is  considered  more  accurate 


than  two-point  interpolation  because  it  does 
not  round  the  origin  point  to  a  gridpoint. 
Eight-point  interpolation  computes  data 
using  every  gridpoint  horizontally  and 
vertically  adjacent  to  it. 

7.16  The  Entrainment  Process.  One 

problem  with  true  Lagrangian  advection  is 
that  there  can  be  no  interaction  between  air 
parcels.  An  entrainment  variable,  therefore, 
has  been  added  to  counter  this  dei^ciency. 
The  model  simulates  two  t3rpes  of  mixing: 

•  The  first  slightly  adjusts  the  CPS  of  the 
parcel  arriving  at  the  termination  point  in 
the  trajectory  using  the  original  CPS  value 
of  that  termination  point.  If  the  original 
CF^  value  is  greater,  the  new  CPS  will  be 
increased;  if  the  original  CPS  value  is  less, 
the  new  CPS  is  decreased.  This  process  is 
applied  to  all  points.  It  tends  to  make  the 
the  difference  between  cloud  and  cloud-free 
areas  less  distinct. 

•  The  second  type  of  mixing  only  applies 
when  ascending  parcels  cause  cloudy  points 
to  become  clear  or  when  descending  parcels 
cause  cloud-free  points  to  become  cloudy. 
The  first  step  is  to  determine  whether  or  not 
the  conditions  causing  these  effects  warrant 
mixing.  If  they  do,  rising  parcels  are  forced 
to  enhance  cloud  formation  and  descending 
paurcels  are  forced  to  inhibit  clouds,  based  on 
(and  proportional  to)  the  surrounding 
conditions  and  the  rate  of  vertical 
displacement. 

7.17  Diurnal  Effects.  5LAYER 
incorporates  diurnal  effects  in  the  cloud 
forecasting  process.  It  assumes  that,  with 
the  right  proportions  of  moisture  and 
heating,  air  parcels  will  rise.  If  they  reach 
the  lifted  condensation  level,  clouds  will 
form.  The  diurnal  effects  variable  takes 
local  effects  into  account,  such  as  those 
caused  by  steep  terrain  features. 
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7.18  Low-L0V«l  Precifirfttttiofi  Evapo* 
ration  »  based  on  two  factors:  (1)  Very  dry 
air  evaporate  at  a  greator  rate  than  moist 
air,  and  (2)  Large  amounts  of  available 
liquid  water  result  in  large  amounts  of 
evapm'ation.  The  amount  of  evapcaration  is 
decreased  whmi  moist  air  is  present  or  when 
available  liq[uid  wat^  is  lacking. 

7.19  Low»LoviH  Molature  Sourets. 
SLAYER  lacks  low*level  moisture  sources, 
such  as  the  Gulf  of  Mexico.  For  this  reason, 

model  has  an  inherent  tendency  toward 
dryness. 


7.^  leing  Forecaat  Variahtoa.  llie  icing 
fcsrecast  product  by  SLAYER  is  based  on 
criteria  in  AWS/TR-SCVOOl,  Forecasters 
Gtdde  on  Aircraft  leing.  Initially,  all  icing 
intensitioB  are  {qpedfied  as  light."  When 
Goid-air  advection  is  forecast  to  ecmeed  2 
hours,  idng  intensity  is  increased  to 
moderate.  When  wartn-air  advection 
exceeds  0.1  K/3  hours,  the  intensity  is 
decreased  to  a  "trace.”  The  criteria  are 
built  into  SLAYER  for  forecasting  icing,  as 
shown  by  Table  4. 


TABLE  4.  SLAYER  Icing  Forecast  Criteria. 
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Chapter  8 

NATIONAL  METEOROLOGICAL  CENTER  PRODUCTION  CYCLES 


8.1  Hie  Aviation  Run  (AVN)  providea 
giol>al  for«casts  of  synoptic  systems  and 
uppe3r-level  winds  out  to  72  hours.  It  is 
performed  every  12  hours,  but  like  the 
AFGWC  GSM  run,  a  6-hour  forecast  is 
produced  as  the  First-Guess  for  the  final 
run,  which  will  be  described  later.  On 
NMC's  CKAY  supercomputer,  the  analysis 
run  takes  about  15  minutes;  the  forecast 
nm,  about  80  minutes.  When  the  CRAY  is 
down,  the  last  forecast  produced,  is  used  as 
a  back-up.  If  the  outage  is  prolonged,  the 
United  Kingdom  model  and  the  AFGWC 
GSM  are  used  as  back-up.  Because  of  time 
constraints,  the  AVN  is  not  run  on  NMC's 
NAS9000  computer.  Transmission  of  the 
AVN  forecast  fields  begins  about  4  hours 
after  database  time.. 

8.1.1  The  Spectral  Statieticeii  interpolation 
Analysis  System  (SSI).  The  SSI  Analysis 
System  is  the  new  analysis  model  used  for 
the  AVN  run;  it  replaces  the  Globa! 
Optimum  Interpolation  (GOl)  analysis 
model  used  previously.  The  analysis  is 
perfoimed  in  the  Global  Data  Assunilation 
System  (GDAS).  This  system  (like  the 
Regional  Data  Assimilation  System,  or 
RDAS)  couples  continuous  data  collection 
with  a  series  of  analyses.  Data  cut-off  time 
is  H+2:45.  The  6-hour  forecast  produced  by 
the  AVN  from  the  previous  model  run  is 
used  as  the  First-Guess.  An  analysis  is 
performed  every  6  hours. 

There  are  some  major  differences  between 
SSI  analysis  and  the  optimum  interpolation 
(01)  techniques  discussed  earlier.  OI 
techniques  analyze  data  by  gridpoint  and  by 
mandatory  level.  Observations  are  estab¬ 
lished  for  each  gridpoint  by  using  conven¬ 


tional  data  (such  as  surface  observations 
and  rawinsonde),  non-  conventional  data 
(satellite  soundings),  and  forecast  data  in  a 
weighted  averaging  system.  SSI,  on  the 
other  hand,  analyzes  ail  data  types  in 
spectral  space  simultaneously.  SSI  is  also 
capable  of  incorporating  Defense 
Meteorological  Satellite  Program  (DMSP) 
Special  Sensor  Microwavq/Imager  (SSMd) 
total  precipitable  water  and  surface  wind- 
speed  data.  The  NMC  Regional  Optimum 
Interpolation  (ROI)  and  the  AFGWC  HIRAS 
analyses  caimot  use  these  types  of  data. 

The  analysis  in  spectral  (as  opposed  to  grid) 
space  hss  two  major  advantages. 

First,  there  is  no  need  to  run  a  normal  mode 
initialization  prior  to  using  the  data  in  the 
forecast  models.  Most  other  operational 
Eissixmlation  systems  require  analyses  to  be 
initialized  because  large-scale  gravity  waves 
are  produced  by  imbeilances  in  the  mass  and 
motion  fields.  lYhen  initialization  is  applied 
to  SSI,  very  small  changes  are  made  in 
contrast  to  the  substantial  alteration  of  Ol 
analyses.  False  \  wes  in  the  tropics,  for 
ex3unple,  are  removed  from  an  OI  analysis 
during  initialization.  The  SSI  analysis  does 
not  produce  those  features,  and  it  is  not 
necesssury  to  attempt  removal. 

The  second  major  advantage  that  the  SSI 
analysis,  accomplished  in  spectral  space, 
uses  levels  close  to  the  operational  forecast 
model  sigma  levels.  In  Ol  analyses,  a  great 
deal  of  interpolation  is  required  to  convert 
the  mandatory  levels  of  the  analysis  to  the 
sigma  levels  of  the  operational  forecast 
model,  resulting  in  analysis  data 
degradation. 


8.1.2  The  NMC  Oiobal  Spectral  Model  (NMC  (this  was  the  GSM  model  acquired  from 

GSM).  The  l^wave  JSMC  GSM  is  the  NMC  by  the  Air  Force).  It  was  increased  to 

forecast  model  for  the  AVN  run.  Resolution  80  waves  in  1988,  and  went  to  the  present 

is  about  equal  to  a  105-km  grid  field  (86  km  126  waves  in  1991.  Its  18  sigma  levels,  as 
less  than  the  AFGWC  half-mesh).  Output  compared  to  the  12  levels  of  the  AFGWC 

fields  are  produced  at  6-hour  intervals  and  GSM,  allows  higher  resolution  in  the  jet  and 

output  on  both  a  65  X  66  point,  381-km  boundary  layers.  An  example  of  the  vertical 

resolution  grid  (the  same  as  the  AFGWC  resolution  of  the  NMC  GSM  is  shown  in 
whole-mesh  grid)  and  in  spectral  coefficient  Figure  21.  The  NMC  GSM  uses  a  mean 
form.  The  GSM  started  as  a  30-wave  model  orography  terrain— see  Figure  22. 
in  1930  and  went  to  40  waves  in  1983 


NMC  MODEL  STRUCTURE 


Flsure  21.  yerticel  Resolution  of  the  NMC  Spectra!  Model  (AVN/MRF),  LFM,  and  NGM  (from 
Hoke  et  ai.,  1989). 


Figtira  22.  Th»  Topography  tho  US,  Canoda  and  Mexico  aa  Recognized  by  the  AVN/MRF 

Model.  The  contour  interval  is  100  meters  (from  the  Technical  Procedures  Bulletin,  1991). 


Mean  orography  terrain  elevations  are  lower 
than  in  the  Nested  Grid  Model  (NGM) 
terrain  (600  to  1,000  meters  lower  over  the 
Rockies,  for  example)  but  it  represents 
reaUty  better  than  the  rolling  hills  of  the 
AFGWC  GSM  terrain.  In  mean  orography, 
elevation  in  the  tropics  is  higher  than  that 
of  the  NGM.  The  I^C  GSM  also  includes 
the  effects  of  gravity  wave  drag,  a  sea- 
surface  temperature  analysis  (the  AFGWC 
GSM  uses  climatological  sea-surface 
temperatures),  dynamically  interactive 
clouds,  and  a  maiine  stratus  parameter¬ 
ization. 


8.2  The  Medium  Range  Run  (MRF). 

The  MRF  is  run  at  OOOOZ  only.  It  provides 
a  medium-range  global  forecast  out  to  10 
days.  Twice  a  month,  it  is  run  out  to  15 
days  on  an  experimental  basis.  The  (XX)0Z 
run  time  was  chosen  to  maximize  the 
amount  of  rawinsonde  data  available  to  the 
model.  The  models  used  for  analysis  and 
forecasting  are  the  same  as  the  AVN  run.  It 
takes  2  hours  and  55  minutes  to  run  the 
entire  analysis  and  forecast  out  to  10  days 
on  the  CRAY.  If  it  fails  on  the  CRAY  after 
the  144-hour  point,  no  backup  is  attempted. 
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If  the  CRAY  is  down  for  a  prolonged  period, 
the  UK  and  AFGWC  GSM  models  are  used 
as  back-up. 

As  with  the  AVN,  MRF  uses  SSI  to  perform 
its  analysis.  A  very  late  data  cut-off  time 
(H-tS  hours)  is  used  to  collect  as  much  global 
data  as  possible.  The  MRF  initialization 
should  be  slightly  better  than  the  AVN's  due 
to  the  increase  in  the  amount  of  data 
available.  Additionally,  late  arriving  data 
firom  locations  such  as  Asia  affecte  the 
forecast  for  the  U.S.  because  of  the  nature  of 
spectral  representation. 

8.3  The  Eaily  Run  (ERL).  The  ERL  is 

run  on  a  12-hour  cycle.  It  provides  a 
quick-look  regional  forecast  update  using  an 
early  data  cut-off,  but  6-  and  12-hourly  data 
continues  to  come  in  long  after  the 
observation  tune.  The  data  cut-off  is  the 
time  observations  must  be  received  in  order 
to  be  used  in  a  particular  run.  The  data 
cut-off  time  for  the  ERL  is  H+l:15. 

8.3.1  The  ERL  Analysis.  The  analysis 
portion  of  the  ERL  run  is  based  on  a  method 
known  as  "successive  correction,"  a  highly 
modified  version  of  the  original  Cressman 
analysis  method  of  1959.  It  is  not  an  01 
anaysis  method.  It  is  produced  by  making 
corrections  to  the  6-hoiir  forecast  fields 
(First-Guess)  generated  in  the  previous  run's 
final  run,  which  will  be  described  later. 
The  corrections  are  interpolated  onto  a 
regularly  spaced  grid  and  added  to  the 
First-Guess  fields  to  produce  the  analysis, 
which  does  not  use  data  from  one  level  to 
modify  data  at  another,  as  an  01  analysis 
would.  Corrections  are  applied  to  the  10 
standard  constant-pressure  levels  through 
100  mb  on  a  129  x  129  point,  190.5-km 
mesh  (equal  to  the  AFGWC  half-mesh)  I^T 
grid. 

8.3.2  The  48-hour  Limited  Fine  Mesh  (LFM). 
The  LFM  is  the  forecast  model  for  the  ERL. 
It  makes  its  forecast  on  seven  sigma  levels: 


a  boundary  level,  three  tropospheric  levels, 
and  three  stratospheric  levels.  The  forecast 
is  run  out  to  48  hours  and  is  the  back-up 
model  for  the  NGM.  There  is  no  back-up 
model  for  the  LFM,  which  is  scheduled  to  be 
turned  off  in  the  near  future. 

When  the  LFM  run  is  complete,  jobs  are 
executed  to  generate  graphic  and 
alphanumeric  bulletins  and  further 
manipulate  the  forecast.  An  example  of  the 
output  is  the  Model  Output  Statistics  (MOS) 
forecast,  which  takes  forecast  elements 
directly  from  the  LFM  outputs  and  uses 
them  as  prediction  variables  in  previously 
derived  regression  equations  to  predict  an 
assortment  of  weather  elements  not  directly 
forecast  by  any  dynamic  prediction  model. 
Examples  include: 

•  Probability  of  precipitation. 

•  Maximum  and  minimum  temperatures. 

•  Indicators  of  severe  weather. 

8.4  The  Regional  Run  (RGL).  RGL  is 
run  on  a  12-hour  cycle  with  a  data  cut-off 
time  of  H'4-2  hours.  It  provides  more 
detailed  regional  forecasts  by  using  full 
rawinsonde  data  and  more  detailed 
precipitation  guidance.  The  RGL  analysis 
and  forecast  ai’e  executed  on  the  CRAY; 
they  each  take  )  6  minutes  to  run.  The 
NAS9()00  is  too  slow  to  run  the  RGL.  When 
the  CRAY  is  down,  the  LFM  is  used  as  a 
back-up. 

8.4.1  Regional  Optimum  Interpolation  (BOI). 
The  ROI  is  used  for  the  RGL  prognosis  (the 
Nested  Grid  Model-NGM).  It  is  similar  to 
the  AFGWC  HIRAS  analysis  in  its  use  of 
optimum  interpolation.  ROI  differs  firom  the 
HIRAS  in  two  ways:  It  does  not  undergo  a 
spectral  transformation  (the  ROI  is  strictly 
a  gridpoint  analysis)  and  it  has  a  different 
run  frequency  than  HIRAS.  ROI  is  run  five 
times  in  a  12-hour  cycle  called  the  "Regional 
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Data  Assiioilation  System,"  or  "RDAS."  T^e 
upper-air  HIRAS  analysis  is  only  run  twice 
in  the  6-hour  AWAPS  cycle.  An  explanation 
of  the  RDAS  cycle  (T-12  to  T-00)  follows. 
The  "T-"  num^rs  represent  data  cut-off 
fimAB  xninus  tihe  number  of  hoinv  indicated 
and  give  the  time  the  analyses  are 


performed.  Analyses  between  T-12  and  T-00 
are  considered  "one-way."  This  means  that 
data  from  outside  the  boimdary  area  can 
affect  the  analysis  within  the  C-grid,  but  not 
vice  versa.  Figure  23  shows  that  portion  of 
the  Northern  Hemisphere  covered  by  the 
C-grid. 


T-12  The  first  analysis  is  performed  using 
oonvenMonal  data  along  with  the  output 
analysis  firom  the  Global  Data  Assimilation 
Ssnitem  (GDAS)  for  the  background  fields; 
this  is  the  same  way  in  which  the  AFGWC 
First-Guess  is  used  in  the  HIBAS  analysis. 
Since  the  prognosis  and  subsequent  analyses 
in  this  cycle  are  primarily  regional,  rather 
than  hemupheric,  initializing  the  boundary 
areas  of  ^e  regional  analyses  by  using  the 
analysis  from  GDAS  is  necessary  to  ensvire 
that  the  longest  wave  components  (which 
extend  beyond  the  regional  grid)  are 
accurate.  When  the  ROI  analysis  is 
complete,  the  NGM  runs  a  S-hoiu*  forecast 
on  the  C-grid  for  use  as  part  of  the 
background  in  the  next  analysis.  This  3- 
hour  forecast  acts  as  a  First-Guess  for  the 
next  analysis.  The  C-grid  is  a  regional  grid 
covering  only  a  portion  of  the  Northern 
Hemisphere. 

T-9  The  second  analysis  is  performed  with 
the  latest  conventional  data,  the  3-hour 
NGM  C-grid  fields,  and  GDAS  data.  The 
data  fiom  GDAS  is  interpolated  from  the 
GDAS  analysis  valid  at  T-12  and  T-6.  This 
interpolated  set  of  fields  is  used  as  boundary 
data  for  the  analysis;  it  only  covers  the  part 
of  the  hemisphere  that  the  NGM  C-grid 
data  does  not  cover.  The  NGM  and  GDAS 
data  is  used  as  the  background  for  the 
analysis. 

T<6  The  third  analysis  is  performed  with 
the  latest  conventional  data,  the  latest 
3-hour  NGM  C-grid  fields,  and  GDAS  data. 
The  data  fix>m  GDAS  is  the  GDAS  analysis 
valid  at  this  time.  The  GDAS  analysis  is 
again  only  used  to  provide  bovmdary 
conditions  outside  the  C-grid  of  the  NGM. 
The  NGM  and  GDAS  data  are  used  as  the 
background  for  the  analysis.  When  the 
analysis  is  complete,  the  NGM  performs 
another  3-ho\ir  forecast. 

T-S  The  fourth  analysis  is  performed 
with  the  latest  conventional  data,  the 


3-hour  NGM  C-grid  fields,  and  GDAS  data. 
The  data  from  GDAS  is  interpolated  from 
the  GDAS  analysis  valid  at  T-S  and  a  6-hour 
forecast  (valid  at  T-OO)  produced  by  the 
AVN.  This  interpolated  set  of  fields  is  used 
as  boundary  data  for  the  analysis;  it  only 
covers  the  part  of  the  hemisphere  that  the 
NGM  C-grid  data  does  not.  The  NGM  and 
GDAS  data  are  used  as  the  background  data 
for  the  analysis.  When  the  analysis  is 
complete,  the  NGM  performs  another  3-hour 
forecast. 

T-00  The  final  ROI  analysis  is 
accomplished  by  using  the  latest 
conventional  data,  the  last  NGM  3-hour 
forecast,  and  the  latest  analysis  data  from 
the  GDAS  to  cover  the  boundary  areas.  To 
accoimt  for  any  differences  between  the 
NGM  3-bour  forecast  smd  the  surrounding 
GDAS  analysis,  an  area  around  (and  ten 
gridpoints  into)  the  C-grid  is  used  as  a 
buffer.  Within  this  area,  the  data  between 
the  NGM  forecast  and  the  GDAS  analysis  is 
interpolated  to  connect  the  two  data  sets 
smoothly.  When  completed,  this  analysis  is 
used  to  initialize  the  run  of  the  NGM's 
48-bour  forecast. 

8.4J2  The  Nested  Grid  Model  (NGM)  is  the 
prognosis  model  used  by  the  RGL.  The  grid 
space  for  the  analysis  is  about  1.5**  latitude 
X  2.0**  longitude  over  North  America,  and 
about  3°  elsewhere.  The  analysis  data 
fields  are  interpolated  to  the  NGM  grids  and 
subjected  to  a  normal  mode  initialization 
that  brings  the  mass  and  motion  fields  back 
into  approximate  balance.  The  vertical 
resolution  for  the  NGM  is  the  same  as  that 
shown  in  Figure  21.  The  forecast,  which  at 
one  time  used  a  nested  set  of  three 
computational  grids  referred  to  as  "A",  "B," 
and  "C,"  now  uses  only  two;  the  A-grid  is  no 
longer  used.  The  B-grid,  with  its  current 
resolution,  has  been  expanded  horizontally 
to  cover  the  entire  Northern  Hemisphere. 
The  C-grid  nas  been  expanded  to  cover  an 
area  a  little  larger  than  the  original  B*grid. 
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Figure  24  shows  the  grids  and  their  areal 
coverage.  The  new  B*grid  provides 
boundary  data  for  the  interior  C-grid,  which 


has  the  finest  horizontal  resolution  of  the 
two. 


Figure  24.  The  NGM  ''Super>C"  Grid  Overlaying  the  NGM  B-Grid.  The  boxes  within  the  C- 
grid  outline  the  areal  coverage  of  the  original  Band  C-gnds  of  the  NGM  prior  to  the  deletion 
of  the  A-grid  (from  Dimego  et  al.,  1991). 
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Tbie  NGM  forecasts  out  to  48  hours,  but  a 
few  elements  (such  as  temperature)  are 
extrapolated  beyond  the  end  of  the  forecast 
period.  The  NGM  forecasts  some  elements 
(sudx  as  vertical  velocities  and  precipitation) 
on  tihe  C’grid.  Hiose  elements  that  do  not 
need  such  a  high  resolution  (500-mb  heights, 
thickness,  surface  isobars,  and  700>mb 
moisturadieights)  are  interpolated  out  to  the 
B-grid. 

Ihe  terrain  used  is  a  mean  orography  with 
greater  peaks  and  valleys  than  l^e  domes  of 
the  AFGWC  GSM  and  the  NMC  Spectral 
Model  (see  Figure  25).  like  the  LFM,  a  hdl 
range  of  MOS  data  is  available  based  on  the 
NGM  forecast.  Since  MOS  data  makes 
corrections  for  a  model's  biases  based  on  its 
historical  performance,  any  signiHcant 
change  to  the  model  intended  to  decrease 
the  effects  of  a  sp>ecific  bias  also  decreases 


the  accuracy  and  reliability  of  the  MOS. 
This  also  applies  to  the  LFM. 

8.5  The  Final  Run  (FNL).  This  run 
produces  the  best  possible  First-Guess  to  be 
used  as  a  foundation  for  the  GDAS  in  the 
subsequent  analysis  cycle.  It  is  run  every  6 
hours.  The  final  run  is  executed  as  a  paired 
analysia/6'hour  forecast  like  the  ATCWC 
First-Guess.  One  difference  is  that  the 
6-hour  forecast  is  the  basis  for  three 
different  analysis  schemes  at  NMC  versus 
just  one  (HIRAS)  at  AFGWC.  The  analysis 
and  forecast  systems  used  are  identical  to 
the  AVh^MRF  runs,  which  have  better 
resolution  than  the  AFGWC  First-Guess. 
The  final  run  takes  20  minutes  on  the 
CRAY.  If  the  CRAY  fails,  the  back-up  is  the 
earlier  AVN  run's  12-hour  forecast.  There  is 
no  attempt  to  perform  intermediate 
(060(yi8(X)Z)  updates. 
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Chapter  9 

STRENGTHS  AND  WEAKNESSES  OF  THE  NMC  MODELS 

Note:  The  meet  important  thing  to  remember  when  considejing  the  inherent  defUAendee  in 
each  of  the  models  is  that  one  problem  can  cause  many  others.  For  example,  a  cold  bias  in  a 
paiiicular  area  usually  causes  heights  to  be  too  low.  This  could  r^ult,  for  just  one  example, 
in  the  jet  b&ng  placed  too  far  south. 


9.1  Limited  Fine  Mesh  (LFM).  The  LFM 
has  been  operational  longer  than  any  of  the 
other  models  discussed  here.  As  a  result, 
information  on  the  LFM's  characteristics  are 
more  accurate  and  complete  than  for  any  of 
the  other  models.  LFM  weaknesses  include 
the  following: 

•  The  LFM  is  too  cold  over  the  Rockies. 

•  The  LFM  has  problems  forecasting 
mrstems  moving  out  of  the  Rockies.  It  often 
lowers  500-mb  heights  and  mean  sea>level 
pressures  too  quickly  east  of  the  Rockies, 
resulting  in  the  mispositioning  of  low-level 
southerly  jets.  Moisture  with  these  jets  and 
associated  surface  fironts  are  moved  too 
quickly  eastward  across  the  plains  and  the 
Mississippi  valley.  The  heaviest  precipi¬ 
tation,  therefore,  is  observed  south  and  west 
of  the  forecast  location. 

•  The  LFM  forecasts  too  many  cyclones 
downstream  from  major  mountain  barriers. 
Arkansas  and  Oklahoma  are  particularly 
vulnerable  to  spuriovis  cyclogenesis. 

•  The  LFM  moves  cyclones  too  slowly  in 
the  fall  and  too  quickly  in  the  spring. 

•  The  LFM  imderforecasts  the  strength  of 
cydones  over  ocean  areas. 

•  The  LFM  rarely  forecasts  cyclogenesis 
adequately  over  Nevada  when  an  upper-air 
trough  digs  towards  Nevada  and  moves  east 
of  W  with  an  upper-level  jet  crossing 
the  Sierra  Nevada  range  between  30  &  40  °N. 


•  The  LFM  is  weak  in  forecasting 
explosive  cyclogenesis  in  the  western  North 
Atlantic.  Ihiring  maximum  deepening,  the 
LFM  was  found  to  lag  behind  by  10  mb  for 
the  first  36-48  hours.  For  all  time  ranges, 
there  is  an  eastward  movement  bias  for 
these  systems.  While  the  12-hour  forecast 
generally  moves  cyclones  too  far 
east-southeast,  the  rest  of  the  prognosis 
forecasts  them  too  far  northeast.  This  also 
causes  cydones  to  move  too  fast,  with  a 
mean  distance  error  of  341  km. 

•  Anticyclones  on  the  48-hour  forecasts  are 
not  strong  enough  (espedally  east  of 
100^  W)  and  are  forecast  too  far  southeast. 
The  error  is  greatest  in  the  eastern  U.S.  and 
on  the  east  coast,  but  it  also  occurs  in  other 
sections  of  the  U.S.  and  off  both  coasts.  The 
mean  distance  error  is  510  km.  Because  of 
the  return  flow  from  these  anticydones, 
advection  of  moisture  and  warm  air  develop 
too  quickly.  This  affects  predpitation 
forecasts. 

•  The  LFM  has  difficulty  forecasting  Arctic 
air  masses.  Sometimes,  when  in  reality  only 
8  low-pressure  trough  forms,  the  sea-level 
pressure  is  lowered  too  much  in  the  cold  air, 
the  1,000-500  mb  layer  is  warmed,  and  the 
surface  low  tracks  too  far  northward  into 
the  cold  air. 

•  The  LFM  has  a  tendency  to  overforecast 
precipitation  in  the  cooler  part  of  the  year 
and  underforecast  it  in.  the  summer. 
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9.2  Th*  Nested  Grid  Moifol  (NGM).  Unrealistic  monaoonal'scale  upalope  flow 

aind  heavy  precipitation  are  often  predicted 

•  The  NCEM  has  a  cold  bias  over  the  Rodby  over  Central  America  for  the  same  reason. 
Motmtains.  The  bias  was  partially  corrected 

by  *  bemiBpbwrir.t«Tnp«i-i>f.iirft  •  Snow  and  ice  fields  are  included  and 

adj\»tm^tOPhillipa  1987).  The  adjustment  updated  weekly  by  tibie  National  Environ- 
iMsened  the  ^ects  of  the  cold  bias  over  the  mental  Satellite  Data  and  Information 
Rooky  Mountains,  but  created  a  warm  bias  Service  (NESDIS).  These  fields  affect  the 
over  the  Appalachians.  The  NGM  is  colder  albedo  in  the  radiation  calculations  as  weU 
than  the  LFM  over  the  Rockies.  as  the  sensible  and  latent  heat  fluxes  at  the 

surface.  For  example,  evaporation  from  a 

•  The  NGM  diurnal  changes  in  frozen  surface  is  not  allowed.  Because  the 

temperature  into  account;  the  LFM  does  ice^anow  fields  used  by  the  NGM  have  a 
not.  coarse  resolution,  enow  cover  aroimd  the 

Great  Lakes  can  be  mistakenly  interpreted 

•  There  is  a  diumal  problem  in  the  NGM's  as  ice  cover,  making  the  lakes  appear  to  be 

ability  to  fcarecast  Twaviirnim  temperatmres  frozen.  This  makes  the  35  gridpoints  over 

in  that  the  surface  maximum  temperature  the  Cheat  Lakes  insignificant  as  moisture  or 

liig«ibftbwd  tbe  iw/wtoVii  irmmlAtiftw  mAmmiiTn  heat  souTces.  The  ability  of  the  mo<fel  to 

by  several  houre.  If  skioi  are  clear  during  forecast  lake-effect  events  decreases  as  snow 

fh«>  day,  the  problem  is  m<HO  pronounced.  accumulates  around  the  lakes  in  early 

As  cloud  cover  increases,  the  problem  winter. 

becomes  less  significant  because  the  model 

recognizes  the  decreasing  effect  of  insolation  *  Since  there  is  orrly  one  snow-cover 
with  increased  cloud  cover.  arralysis  during  the  week  (usually  on 

Morrday),  significant  changes  in  mow  cover 

•  There  is  a  pronounced  diurnal  problem  during  the  week  will  not  be  noted.  Since 

with  low-level  wind  forecasta  in  that  the  NGM  caruiot  melt  anow,  the  surface 

fmecasts  in  the  lowest  35  mb  (the  midpoint  temperature  cannot  be  above  0°  C  in  an 

is  shout  160  meters)  peak  at  night,  with  a  area  of  thin  snow  cover,  even  when  there  is 

minimum  at  midday.  In  order  to  produce  a  strong  insolation.  The  existence  of  snow 

more  realistic  forecast  at  anemometer  level,  cover  in  the  model  at  a  point  where  it 

the  model  needs  higher  verticsd  resolution  should  be  mow-free  hsiB  been  shown  to 

f>nd  a  more  realistic  method  for  modeling  reduce  the  bottom-layer  sdr  temperature  in 

the  downward  transport  of  momentiun  some  cases  by  over  6^  C.  When  forecsutting 

during  the  day.  Currently,  the  model  rain  or  snow,  therefore,  forecasters  must 

reflects  a  wind  more  typicsd  of  a  nocturnal  consider  the  model's  deficiencies.  Note  also 

jet.  that  heights,  especiaTy  in  the  mid-levels, 

could  be  forecast  too  low  or  too  high 

•  Due  to  its  horizontal  resolution  (86  km),  depending  on  the  mow  cover,  since  the 

NGM  tends  to  forecast  sea  breezes  several  albedo  or  the  effects  of  insolation  may  be  the 

hundred  kilometers  farther  inland  than  they  opposite  of  what's  really  happening, 

actually  are.  Due  to  the  strong  convergence 

forecast  so  far  inland,  the  NGM  also  *  The  NGM  has  a  tendency  to  forecast  the 

overforecasts  the  extent  of  sea-breeze  movement  of  sirrface  cyclones  too  slowly.  In 

enhanced  precipitation  in  those  areas.  The  a  study  by  Grumm  and  Siebers  (1989),  the 

effect  also  occurs  around  the  Great  Lakes.  largest  positional  errors  were  found  in  the 
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southeast  quadrant  of  the  NGM.  Forecast 
positions  w«re  found  to  be  northwest  of  the 
analyzed  position.  The  mean  distance 
errors  grew  from.  90  km  at  the  24-hoiir  point 
to  310  km  at  the  48*hour  point.  The  second 
largest  error  was  found  in  the  southwest 
quadrant.  Cyclones  were  too  far  nortii  of 
the  analyzed  position.  The  eiror  increased 
&om  about  160  km  at  the  24-hour  point  to 
200  km  at  the  48-hour  point.  No  bias  was 
detenoined  in  the  northern  quadrants, 
because  the  forecast  versus  analyzed 
positions  were  too  scattered. 

•  Grumm  and  Siebera  determined  that  the 
NGM  underdevelops  lows  in  the  Atlantic  off 
the  New  England  coast;  this  is  especially 
true  of  1}omb8"  that  begin  development  near 
Cape  Hatteras.  However,  NGM  low 
predictions  are  too  deep  over  the  North 
American  continent  at  all  times.  A  checklist 
was  developed  to  forecast  "bombs"  off  the 
east  coast.  In  most  cases  a  vortidty 
maximum  moves  over  the  already  strong 
thermsd  gradient  between  the  Labrador 
current  and  the  Gulf  stream.  Cyclogenesis 
occurs,  forr  .ing  a  new  low  or  deepening  am 
already  e*  ^sting  one.  The  low  then  deepeiui 
dramatically  as  it  moves  N  and  NE.  The 
checklist  developed  by  Sanders  and  Auciello 
(1989)  applies  only  to  the  Atlantic  north  of 
40®  N  and  west  of  60®  W.  The  checklist 
predicted  20  "bombs"  between  October  1^8 
end  March  1989  and  correctly  forecaist  13  of 
the  16  that  occurred.  At  least  four  of  the 
following  six  criteria  must  be  met; 

1.  The  strength  of  the  500-mb  absolute 
vorticity  maximum  is  equal  to  or  greater 
than  1.7  X  s  '. 

2.  The  vorticity  maximum  stays  the  same  or 
strengthens  during  the  NGM  forecast. 

3.  The  vorticity  maximum  is  forecast  to 
move  faster  than  30  knots. 


4.  The  vorticity  maximum  crosses  the 
Atlantic  coast  between  32  and  45“  N. 

5.  There  is  a  jet  max  greater  than  100 
knots  at  250  mb  or  3(X>  mb,  just  south 
of  tiie  initial  vorticity  maximvun. 

6.  The  NGM  forecasts  a  suHace  low  of 
990  mb  or  deeper  over  the  Atlantic 
Ocean  north  of  38“  and  west  of  55® 

•  The  NGM  often  over-deepens  lows 
coming  out  of  the  Rocky  Mountains. 

•  Another  study  by  Grumm  suid  Siebera 
(for  the  period  November  1988-January 
1989)  found  that  anticyclones  were 
incorrectly  positioned  by  261  km  and  334 
km  at  the  24*  and  36-hoiir  points, 
respectively.  The  latter,  however,  is  much 
better  than  the  LFM's  510  km  error.  The 
NGM  also  slightly  over-forecasts  pressure 
centen  in  the  Rockies  and  undm^fo^ecasts 
them  in  the  eastern  U.S.,  probably  due  to 
the  cold  and  warm  biases,  respectively. 

•  NGM  precipitation  forecast  amounts  are 
superior  to  the  LFM  and  NMC  GSM  model 
when  the  associated  500-mb  trough  is 
positively  tilted. 

•  The  NGM  has  a  problem  forecasting 
heavy  amounts  of  precipitation  (greater 
than  .6  inch/day)  in  the  southern  U.S. 
during  the  cool  season.  A  study  by  N.  W. 
Junker  (for  the  winters  of  1987-88  and 
1988-89)  showed  the  greatest  number  of 
forecasting  misses  to  be  along  the  GxUf 
coast,  but  decreasing  northward. 

•  The  NGM  over-forecasts  light  precipi¬ 
tation  in  the  summer. 
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9^  Tli«  NMC  GSM  Model.  Most  of  the 
error  discussions  that  follow  should  be 
interpreted  with  caution,*  they  are  averaged 
over  a  limited  period  of  tame  and  may 
depend  on  the  flow  regime.  Also,  the  GSM 
was  recently  upgraded  from  an  SO*  to  a 
126*wave  modbl. 

•  Ihe  80-wave  GSM  had  a  general  (»ld 
bias  that  was  greatly  reduced  by  the  new 
126-wave  model,  which  appears  to  have  a 
slight  cold  bias  in  the  mid-latitudes. 

•  The  new  1^-wave  model  appears  to 
handle  coastal  and  lee-side  cyclogenesis 
better  dum  the  80-wave  version  and  better 
than  the  NGM  and  LFM  models. 

•  Ihe  GSM  model  has  an  upper-level 
easterly  bias  over  the  tropics. 

•  In  summer,  heights  are  over-forecast 
over  land  and  under-forecast  over  oceans. 
These  erton  have  been  linked  to  apparent 
thermal  forcing. 

•  The  GSM  under-forecasts  the  intensity 
of  Arctic  highs  in  centi*al  and  western 
Canada. 


9.4  The  Spectral  Statistical  Inter¬ 
polation  Analysis  System.  Due  to  the 

relative  xmwness  of  this  model,  information 
on  its  peiformance  is  limited.  Testing, 
which  concentrated  on  the  system's  effects 
on  forecast  models,  showed  slight  overall 
improvement  in  forecast  skUl  up  to  day  5 
when  compared  to  a  parallel  forecast 
initialized  by  the  old  OI  analysis.  Beyond 
day  5,  no  discernible  differences  were  noted. 


•  The  GSM  makes  the  plunge  of  Arctic  air 
into  the  U.S.  too  slow. 
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GLOSSARY 

. 

APGWC 

Air  Global  Weather  Central 

-  ■ 

» 

ASDAE 

Aircraft  to  Satellite  Data  Relay,  a  syst^  m  which  wide^bodied  aircraft  are 
fitted  with  equipment  to  obtain  meteorological  data. 

AVN 

Aviation  run.  National  Weather  Servi<» 

AWAPS 

Advanced  Weather  Analysis  and  IVediction  System,  Air  Force  Global  Weather 

Central 

AWSf^ 

Air  Weather  Service  Primitive  Equation,  the  forecast  model  used  at  Air  Force 

Global  Weather  Central  prior  to  acquisition  of  the  Global  Spectral  Model  by  the 

Air  Force 

CPS 

Condensation  Pressure  Spread 

CPU 

Computer  or  Central  IVocessing  Unit 

DMSP 

Defense  Meteorological  Satellite  Program 

ERL 

Early  Run,  National  Weather  Servi<» 

• 

FNL 

Final  Run,  National  Weather  Service 

GADB 

Global  Applications  Database,  Air  Force  Global  Weather  Central 

GDAS 

Global  Data  Assimilation  System  *  National  Weather  Service 

GOI 

Global  Optimum  Interpolation,  National  Weather  Service 

GSM 

Global  Spectral  Model,  Air  Force  Global  Weather  CentrsJ  and  National  Weather 

Service 

HIRAS 

High  Resolution  Analysis  System,  Air  Force  Global  Weather  Central 

LFM 

limited  Fine  Mesh,  National  Weather  Service 

MRP 

Medium  Range  Forecast  run.  National  Weather  Service 

NGM 

Nested  Grid  Model,  National  Weather  Service 

' 

NESDIS 

National  Environmental  Satellite,  Data,  and  Information  Service 

NMI 

Normal  Mode  Initiedization 

• 
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NWS 


National  Weather  Service 


01  Optimiim  Interpolation 

Polar-stereograghic 
QC  Quality  Control 

RAOB  Radi(M!ond^Rawinsonde  Observation 

HDAS  Begional  Data  Assimilation  System,  National  Weather  Service 

RGL  Regional  Run,  National  Weather  Service 

ROl  Regional  Optimum  Interpolation,  National  Weather  Service 

RTNEJFH  Real'Time  Nephanalysis  model,  Air  Force  Global  Weather  Central 

RWM  Relocatable  Window  Model,  Air  Force  Global  Weather  Central 

SDKS  Satellite  Data  Handling  System,  Air  Force  Global  Weather  Central 

SGDB  Satellite  Global  Database,  Air  Force  Global  Weather  Central 

SNODEP  Snow  Depth  model,  Air  Force  Global  Weather  Central 

SSI  Spectral  Statistical  Interpolation,  National  Weather  Service 

SSj^VI  Special  Sensor  Microwave/Imager,  a  Defense  Meteorological  Satellite  Program 
sensor 

SLAYER  The  primary  cloud  forecast  model  used  at  Air  Force  Global  Weather  Central 
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APPENDIX  A 

APGWC  SATELLITE  DATA  HANDLING  SYSTEM  <SDHS)  PRODUCTS 


AFGWC'produced  SDH5  facsimile  products  are  listed  below,  along  with  the  models  used  to 
support  them. 


Facsimile  Chart 

1.  ASNH 

2.  FSNH2 

3.  FSNH3 

4.  FUNH2 

5.  FUNH3 

6.  FHNH2 

7.  FHNH3 

8.  FANH25 

9.  FANH3 

10.  WWUS12 

11.  WWUSA-C 

12.  FANA 

13.  FENH66-69 

14.  FENHIA 

15.  FATR 


Numeric  Model 

First-Guess,  HIRAS,  5LAYER. 

AFGWC  GSM,  supplemented  with  NMC  GSM  (AVN  RUN)  and  NGM. 

AFGWC  GSM,  supplemented  with  NMC  GSM  (AVN  RUN)  and  NGM. 

6LAYER.  AFGWC  GSM,  supplemented  with  NMC  GSM  (AVN  RUN) 
and  NGM. 

SLAYER,  AFGWC  GSM,  supplemented  with  NMC  GSM  (AVN  RUN) 
and  NGM. 

SLAYER,  AFGWC  GSM,  supplemented  with  NMC  GSM  (AVN  RUN) 
and  NGM. 

SLAYER,  AFGWC  GSM,  supplemented  with  NMC  GSM  (AVN  RUN) 
and  NGM. 

LAYER,  AFGWC  GSM,  supplemented  with  NMC  GSM  (AVN  RUN) 
and  NGM. 

SLAYER,  AFGWC  GSM,  supplemented  with  NMC  GSM  (AVN  RUN) 
and  NGM. 

AFGWC  GSM,  supplemented  with  NMC  GSM  (AVN  RUN)  and  NGM. 

AFGWC  GSM,  supplemented  with  NMC  GSM  (AVN  RUN)  and  NGM. 

SLAYER.  AFGWC  GSM,  supplemented  with  NMC  GSM  (AVN  RUN) 
and  NGM. 

NMC  GSM  (MRF  RUN) 

NMC  GSM  (MRF  RUN) 

AFGWC  GSM 
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16.  FAME 


AFGWC  GSM 


APPENDIX  B 


AFGWC  COMPUTER-GENERATED  PRODUCTS 


F*c«imii0  Chart 

Numeric  Model 

Faceimile  Chart 

Numeric  Model 

1.  AXASA-DKGWC 

HIRAS 

12.  FUEUKGWC 

AFGWC  GSM 

2.  AXIO(j)KGWC 

KIRAS 

13.  FUNH(u)KGWC 

AFGWC  GSM 

3.  AXPWKGWC 

HIRAS 

14.  FUNT(ii)KGWC 

AFGWC  GSM 

4.  FAEHKGWC 

AFGWC  GSM 

15.  PUPA(ii)KGWC 

AFGWC  GSM 

5.  FAEH18KGWC 

AFGWC  GSM 

16.  FUTR(ii)KGWC 

AFGWC  GSM 

6.  FAWHKGWC 

AFGWC  GSM 

17.  FUWH(ii)KGWC 

AFGWC  GSM 

7.  FEAS50KGWC 

AFGWC  GSM 

18.  FUXP(ii)KGWC 

AFGWC  GSM 

8.  FEEUKGWC 

AFGWC  GSM 

19.  FXAS(i)KGWC 

AFGWC  GSM 

9.  FENHKGWC 

AFGWC  GSM 

20.  FXEU(i)KGWC 

AFGWC  GSM 

10.  FSNTKGWC 

AFGWC  GSM 

21.  FXICXDKGWC 

AFGWC  GSM 

11.  FUCAKGWC 

AFGWC  GSM 

22.  FXPW(i)KGWC 

AFGWC  GSM 
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APPENDIX  C 

NMC  COMPUTER-GENERATED  PRODUCTS 

NMC-produced  facsimile  products  are  listed  below,  along  with  the  models  used  to  support 
them. 


Facsimile  Chart 


1.  ANAL  600  mb  WITH  PLOT  DATA 

2.  ANAL  700  mb  WITH  PLOT  DATA 

3.  ANAL  860  mb  WITH  PLOT  DATA 

4.  OBS  WND  ALOFT  4-PANEL 

6.  ANAL  700  mb  HGT/MEAN  RH  SFC-aprox  460  mb 

6.  ANAL  850  mb  HGT/ISOTHERM 

7.  ANAL  600  mb  HGT/REL  VORT 

8.  ANAL  MEAN  SLP/lOOO-eOO  THKNS 

9.  ANAL  200  mb  TROPIC  WITH  DATA 

10.  COMPOSITE  MOISTURE  CHART/FREEZING  LVL 

11.  COMPOSITE  MOISTURE  CHAKT/SFC-SOO  mb  MEAN  RH 

12.  COMPOSITE  MOISTURE  CHARTAdFTED  AND  K  INDICES 


Numeric  Model 


EEL  ANAL 
ERL  ANAL 


EEL  ANAL 


LFM 


ROI  ANAL 
ROI  ANAL 


ROI  ANAL 


ROI  ANAL 
ERL  ANAL 
ROI  ANAL 
ROI  ANAL 


ROI  ANAL 


13.  COMPOSITE  MOISTURE  CHART/SFC.600  PRECIPITABLE  WATER  ROI  ANAL 


14.  ANAL  300  mb  WITH  DATA 


16.  ANAL  SFCyiOOO-600  mb  THKNS 

16.  ANAL  600  mb 

17.  ANAL  600  mb  MERCATOR  PROJECTION 

18.  ANAL  700  mb  MERCAIDR  PROJECTION 

19.  PROG  36  hr  300  mb  HGT/WND 


ERL  ANAL 


SSI  ANAL 


SSI  ANAL 


SSI  ANAL 


SSI  ANAL 


NMCGSM(AVN  RUN) 
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20.  PROG  36  hr  700  mb  HGT/WND 

21.  PROG  36  hr  700  mb  HGWND 
JH.  PROG  86  hr  500  mb  HGT/WND 

23.  PROG  12  hr  WND  ALOFT  8  PANEL 

24.  PROG  12  hr  700  mb  HGT/MEAN  RH  SFC-aprox  450  mb 

26.  PROG  12  hr  850  mb  HGT/ISOTHBEM 

26.  PROG  12  hr  500  mb  HGT/REL  VOKT 

27.  PROG  12  hr  MEAN  SLP/1000-600  THKNS 

28.  PROG  24  hr  700  mb  HGIVMEAN  RH  SFC-aprox  450  mb 

29.  PROG  24  hr  850  mb  HGT/Ii^OTHERM 

30.  PROG  24  hr  500  mb  HGT/REL  VORT 

31.  PROG24hrMEANSLPA000-600THKNS 

32.  PROG  36  hr  700  mb  HGT/MEAN  RH  SFC-aprox  460  mb 

33.  PROG  36  hr  850  mb  HGT/ISOTHERM 

34.  PROG  36  hr  500  mb  HGT/REL  VORT 

35.  PROG36hrMEANSLPA000-600THKNS 

36.  PROG  48  hr  700  mb  HGl/MEAN  RH  SFC-aprox  450  mb 

37.  PROG  48  hr  860  mb  HGT/ISOTHERM 

38.  PROG  48  hr  600  mb  HGT/REL  VORT 

39.  PROG  48  hr  MEAN  SLPAOOO-500  THKNS 

40.  PROG  48  hr  500  mb  HGT/WND 

41.  PROG  36  hr  500  mb  HGT/WND 

42.  PROG  2-PANEL  24  &  36  hr  SFC/THKNS 

43.  PROG  4-PANEL  ANAL,  12,  24  &  36  hr  (BAROCUNIC) 


NMC  GSM  (AVN 
NMC  GSM  (AVN 
NMC  GSM  (AVN 


NMC  GSM  (AVN 
NMC  GSM  (AVN 
NMC  GSM  (AVN 
NMC  GSM  (AVN 


RUN) 

RUN) 

RUN) 

LFM 

NGM 

NGM 

NGM 

NGM 

NGM 

NGM 

NGM 

NGM 

NGM 

NGM 

NGM 

NGM 

NGM 

NGM 

NGM 

NGM 

RUN) 

RUN) 

RUN) 

RUN) 
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44.  PROG36hrSFCyrHKNS 

46.  PROG  24  hr  SFC/THKNS 

46.  PROG  4-PANEL  PRECIPHlOBABttJTY 

47.  HiOG36hr700mb 

48.  PROG  48  hr  low  mb  MERCATOR  PROJECTION 

49.  PROG  36  hr  700  mb  MERCATOR  PROJECTION 

50.  PROG  36  hr  BOUNDARY  LAYER  WINDER 

51.  PROG  4-PANEL  24,36.48  &  60  hr  MA15/MIN  TEMPS 

52.  PROG  84  hr  500  mb  HGTS 


NMCGSM(AVN  RUN) 
NMCGSM(AVN  RUN) 
LFM 

NMCGSM(AVN  RUN) 
NMCGSM(AVN  RUN) 
NMCGSM(AVN  RUN) 
NMCGSM(AVN  RUN) 
NGM 

NMCGSM(MRF  RUN) 
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9  OSETOW,  Gruman  Ave.,  Bldf  1060,  Beale  APB,  CA  96903-5000  . I 

96  OSSmW.  Baae  Ope  Rd.,  Bldf  9001,  Dyeaa  AFB,  TX  790074000  . I 

5  OSEDOW,  221  Filfht  Line  Dr..  Bldf  746,  Minot  AFB.  MD  58705-5021  . .  1 
384  OS^DOW,  Kanaaa  Ct,  Ste  104,  Bldf  1112,  McConnell  AFB, 

KS  67221-5000  . I 

2  WECC,  Bldf  6646,  245  Davie  Ave,  Barkadale  AFB,  LA  7111(k500t . 1 

42  OSEDOW,  Oeorfla  Rd..  Bldf  8200,  Rm  10,  Lortnf  AFB.  MS  047514000  1 
416  OSEOSW,  582  Batifar  Rd.,  Ste  121,  Bldf  1000,  Grifliaa  AFB, 

NY  134414000  . 1 

380  OSEDOW,  Alabama  Ave.,  Bldf  2712,  Rm  100,  PUttaburfh  AFBt 
NY  129036000  . I 

97  OSEDOW,  Baker  AFB.  AH  72317-6000  . 1 

319  OSEDOW,  865  Steen  Ave.,  Bldf  528,  Ste  106,  Grand  Fbrka  AFB. 

ND  582064244  . 1 

28  OSECSAW,  Bldf  7606,  Bllaworth  AFB,  SD  577064000 . ) 

361  OSEDOW,  Bldf  S3S,  Whiteman  AFB,  MO  853054000  . 1 

7  OSEDOW.  Bldf  142^  Carewall  AFB,  TX  761274000  . 1 

410  OSEDOW,  F  Ave.,  DMf  401,  Ste  7,  K.l.  Sawyer  AFB.  Ml  49M33400  - 1 

305  OSEDOW,  Heoaler  Blvd.  BMf  S-9B.  Qriaann  AFB,  IN  468714000 - I 

379  OSEDOW,  Oparatlana  St.  BMf  PI4,  Wurtarallh  AFB,  Ml  487534000 . 1 
ACODOW.BIdf  2I,^0BlmStStatl^UntlayAFB.VA2^853««3  ...I 


DIST-1 


WWS,  Unit  0840,  APO  AA  340O1-S0O0 . 1 

MXi^AOSW,  BU(  1130,  Shaw  Dr.,  Shaw  AFB,  SC  29I62-S410 . ] 

QUA.  aeSPW/DOMS,  RFD  1,  McEntin  ANGB,  Baalnvar,  SC  29044-0600  . .  1 
4  Oparationa  GroupOSW,  1080  Curttaa  Ave,,  Sta  140,  Bldf  4607, 

Saymour  Johaaon  APB,  NC  27631-2624  . 1 

31  OSaOSW,  Coral  Saa  Bird.  Bid#  702,  Homaataad  AFB,  FL  33830-6010  .  I 

1  OCVWS,  Fll«htllna  Ara..  Bld(  764,  Lantlay  APB,  VA  23006-6606 . 1 

347  OSSG6W.  Sauannah  St.,  Bld«  822,  Moody  APB,  GA  31000-6000  . I 

68  OSaOSW,  Bld«  P3,  Rn  101,  Florida  Ava.,  McDill  APB,  PL  33808-6000 . 1 

12AO(VAOSW.  B  Ara.,  Bld«  2000,  Barfitnm  APB,  TX  78743-6000  . I 

27  08V0SW,  110  B  Sextant  Ava.,  Sta  1040,  Cannon  APB,  NM  88103-6322  1 
366  OSMMWP,  Phoanix  St,  Bld«  4820,  DarU-Monthan  APB, 

AZ  86707-8801  . 1 

40  os^oew,  BMc  671,  Holloman  APB,  NM  88330-6000  . 1 

68  Oe«OSW,B(h  St.  BMi  063,  Luka  APB,  AZ  86300-6000  . 1 

OL-K  680SaOeSW,  Bide  324,  GiU  Band  APAP,  AZ  85337-6000  . I 

67  06BOSW,  27  Depot  Rd.,  Bide  M6,  Nallia  APB.  NV  80101-6000 . I 

388  0Sa<O8W,  866  Oak  St,  Mt  Homo  APB,  ID  83848-6000  . 1 

24WB<CC.  APO  AA  34001-6000 . 1 

ATSVCDW,  US  Amy  Inul,  Ft  Huaehuea  AI,  AZ  868136000 . 1 

120  ARRSIXIV,  NAS  MoAeU  Fid.  CA  040364008  . 1 

USCENTCOHCCJ3W,  MaoDill  AFB.  PL  338037001  . ) 

USSOOCEN1/SCJ2-SWO.  MaoDill  APB,  PL  338037001 . 1 

USSOCOMBOJ3W.  MacDUl  AFB,  FL  338038001  . I 

Dat  1,  NEADa'DOW,  Mainaiao  Ava.,  Bldf  610,  Banfor  ANGB,  ME 

04401-3000  . 1 

101  AREFW/DO,  Banrw  ANGB,  ME  04401-3009  . 1 

IWG,  Bld«  188,  Haidoa  St,  Ft  MePharaon,  GA  303036000  . 1 

OL-A.  6WS,  APZJ^WO,  Ft  Irwin,  CA  023133000  . ) 

Dat  1.  IWXG,  Bide  7183,  Ft  Campball,  KY  42223-6000 . 1 

Dat  2,  6WS.  Daviaon  USA  AlrAald,  Ft  Balrolr,  VA  220836748 . 1 

OL-B,  Dat  2,  6WS,  Ft  Indiantown  GAP.  PA  170038002  . 1 

Dat  3.  IWXG,  Pra«ar  St,  Bldf  AT3661.  Ft  Bra«f  AIN,  NC  28307-6000  . . .  1 

CDR  USASOC,  Attn:  AOIN-ST.  Ft  Brace.  NC  283078200  . 1 

JSOCVWeaUiar.  P.O.  Box  70230,  Ft  Brace.  NC  28307-6000 . 1 

Dat  VWXG,  Haiiear  Amaaa  Dr.,  Bide  2006.  Ft  Drum,  NY  13802-6042  - 1 

Dat  6,  IWXG.  Bide  6220.  PO  Bax  040,  Pt  Knox.  KY  40121-6640 . 1 

Oat  8,  IWXG,  Bide  3082,  AirpoK  Dr.,  Pi  Lewie,  WA  98433-6000  . 1 

Dat  7,  IWXG,  Bide  618,  Ritehia  AAF,  Pt  Ord,  CA  93041-6111  . 1 

Dat  8,  IWXG,  743  Ray  Plane,  Bide  743,  Pt  Rilay,  K8  884428317  . 

OL-A.  Dat  8AVXG,  Foumay  AAF,  Bide  M04,  Pt  Leonard  Wood, 

MO  86473-6882  . 

Dal  9.  6W9.  Ft  Rucker,  AL  363628182  . 

Det  10,  IWG,  BMc  2486.  Rm  110,  Lamon  AAF,  Ft  Banninc,  GA  310068034 

7Sch  RCR  (Attn:  SWO),  Ft  Banninc,  GA  31906-6000 . 

Det  11.  IWXP,  Pool  Rd.,  Bide  4907.  Ft  Sill,  OK  73503-5100  . 

Det  12,  IWXG,  Bide  880,  Pt  Darana,  MA  01433-6310  . 

Dat  13,  IWXG,  Condon  Rd.,  Bide  2408,  Pt  Euatia,  VA  238048252  . 

Det  14.  IWXG,  BMe  90049.  Rm  1 18,  Gray  AAF,  W  Ft  Hood.  TX  76644-6078 
OLA,  Oat  14.  IWXF,  Slawitxke  St.  Bide  11210.  Ft  Bllaa,  TX  79016-2418  . 
OLA,  Oat  21. 1  WXC,  Straneh  St  Bide  1262.  Hunter  AAF,  CA  31409-6193 

Dat  31.  IWXG,  Pt  Wk  AFF,  Bide  E928,  Ft  Pblk,  LA  714608280  . 

Det  68,  IWXG,  BMe  B601,  BuUa  AAF,  Ft  Caiaon,  CO  800138403  . 

H4  5lh  U.S.  Amy,  AFKB-OP  (SWO),  Ft  Sam  Houeton,  TX  78234-7001  . . 

USAFE(DOW.  Unit  3060,  Box  16.  APO  AE  09094-5000 . 

USAFBTOWO,  Unit  3060,  Box  600,  APO  AE  00094-6000  . 

17AF/DOW.  Unit  4086,  APO  AE  00138-5000 . 

Det  1,86  OPS  OP,  Unit  8496,  APO  AE  00004-6000  . 

OLC,  Dal  1,88  OPS  GP,  Unit  3665,  APO  AE  00122-5000  . 

OLD,  Det  1,  80  OPS  CP,  Unit  7890,  APO  AE  09122-6000  . 

HQ  USEUCOM  ECJ3.  Unit  30400,  Box  1000,  APO  AE  09128-4209  . 

100  OSSiDOW,  Unit  4946,  APO  AE  00459  . 

48  OSSWX,  Unit  8300,  APO  AE  09484-6000  . 

36  TPW,  Unit  3896,  APO  AE  00132-6000  . 

65  AL/8S,  APO  AE  00720-6000  . 

7WS,  CINCU8AREUIVAREAWX,  APO  AE  00403-6000 . 

7WS,  Unit  29361,  APO  AE  09014-6000  . 

OLC,  7WS,  Cnu-  446,  Box  280,  APO  AB  00046 . 

OLD,  7WS.  Unit  23831,  APO  AE  00180-2784  . 

OLE,  7WS,  APO  AE  00137-6707  . 

OLF,  7W8.  Unit  31401,  Box  6.  APO  AE  09830-6000  . 

OLO,  7W8,  Unit  23848,  APO  AE  00272-6000  . 

OLI,  7W8,  Unit  38703,  APO  AE  00236-6000  . 

OLJ.  7WS,  CMR  431,  APO  AE  00175-6000  . 

Det  1,  7W8,  Unit  27210,  APO  AE  00002-0216 . 

Dal  2, 7W8.  Unit  2300,  APO  AE  00166-0816 . 

Dat  3. 7W8,  Unit  20231,  APO  AE  00102-3737  . 

Del  4, 7W9,  Unit  31030,  APO  AE  000264)261 . 


Det  6.  7WS.  Cmr  463,  APO  AE  00146-0079 . 1 

Det  7,  7WS,  Unit  28130,  APO  AE  00114-6000 . 1 

OLA.  Det  7,  7W8,  Unit  28216,  APO  AB  00173-6000 . 1 

Det  8,  7WS,  Unit  26302,  APO  AB  00070-6000 . 1 

Det  9,  7W8,  Cmr  423,  APO  AB  00107-0200 . 1 

OLB,  7WS,  Cmr  422,  APO  AE  00274-6000  . I 

OLA.  Dat  10,  7W8,  CMR464,  Unit  31020,  APO  AE  00260-  0047  . 1 

Dat  10,  7WS,  Unit  28410,  APO  AE  001820006 . I 

OLB.  Dat  10,  TWS,  APO  AE  00031-6000  . 1 

Det  11,  TWS,  Coleman  Bka  Waa,  Unit  29710,  APO  AE  00028-3728 . 1 

Det  12,  TWS,  Unit  24220,  APO  AE  00186-6000 . 1 

OLA,  Det  12,  TWS,  APO  AE  00111-6000  . 1 

Dat  13,  TWS,  Cmr  416,  APO  AE  00140-9008 . 1 

Det  26,  TWS,  Unit  20632,  APO  AE  00006-5000 . I 

4360S5VD0W,  APO  AE  00097-6000  . 1 

81  TFW/WE,  Unit  6075,  APO  AE  00497-6000  . 1 

20OSaiDOM,  Unit  6476,  APO  AE  09468-5000 . 1 

10  TFW/DOM,  Unit  6686,  APO  AE  00470-6000 . I 

320S^WE,  Unit  6706,  APO  AE  00719-6000  . I 

8eWF,  Unit  3000,  APO  AE  00004-6000 . 1 

40  SW/WX,  Unit  6160,  APO  AE  00801-6000  . 1 

62  OS^WEP,  Unit  3720,  Bax  196,  APO  AE  00126-6000  . 1 

16AF/SWO,  Unit  6366,  APO  AE  00841-6000 . 1 

401  CSCVWE,  Unit  8400,  APO  AE  00641-6000  . 1 

39  TACeVDOM,  Unit  7000,  Box  1  IS,  APO  AE  00824-6000  . 1 

62TFW/OTM.  APO  AE  00706-5000  . 1 

OLC.  TWS,  CMR445,  Box  280,  APO  AE  09046-5000  . 1 

48  Weather  Flieht,  APO  AE  00464-6000  . 1 

106  Weather  Flieht,  Tenneaeaa  A>  National  Guard,  PO  Box  17267, 

Naahvilla.  TN  37167-2001  . 1 

107  Weather  Flieht,  Seirridee  ANGB,  Ml  48046-6024  . 1 

1 10  Weather  Flieht,  MO  ANG,  131TFW,  Brideeton,  MO  03044-2371  . 1 

111  Weather  Flieht,  Bllincton  ANGB,  TX  77034-5686 . 1 

113  Weather  Flieht,  Hulman  Fid,  Tatra  Haute,  IN  47830-6000 . 1 

116  Weather  Flight,  WA  ANG.  Bldg  304,  MeChord  AFB,  WA  98433-6000  . .  1 

120  Weather  Flieht.  Buckley  ANGB.  CO  80011-0600  . 1 

121  Weather  Flieht.  8to|i  28,  Andtm  AFB,  MD  30331-8630  . 1 

122  Weather  Flieht.  New  Oriaana  NAS,  LA  70143-0200  . 1 

123  Weather  Flieht,  Pnrtland  lAP,  OR  07218-2707  . 1 

125  Weather  night,  PO  Box  680340,  Tulaa  AFS,  OK  741684)340 . 1 

128  Waathar  night,  WIANC,  360  E  Collaca,  Milwukae,  Wl  53307-0208  ...  1 
127  Waathar  night,  Forbaa  nd.  Topaak,  KS  8e619«)00  . 1 

130  Waathar  night,  Yeager  Apt,  Chariaaten,  WV  2631 1-6000 . 1 

131  Weather  night,  Bamaa  Map.  WaatOaM,  MA  01085-1386  . 1 

140  Weatlwr  night.  Willow  Gruva  NAS,  PA  19000-6106 . 1 

146  Weather  night,  CTR  Plllaburg  ANG  AN,  PA  162314)460  . 1 

154  Waathar  night.  Camp  Rnbinaon,  North  Little  Rock,  AR  72118-2300  . .  1 

168  Waathar  night,  6226  Morria  nd  Dr.,  (nwrlotle,  NC  28208-5707  . 1 

160  Waathar  night,  efo  HQ  FLANG,  Suia  Araanal,  St  Auguatina, 

FL  32086-1008  . 1 

164  Weather  night,  Rkkenbaekar  ANGB,  OH  43217-6007  . 1 

166  Weather  night,  SUndiford  FM,  Loulavilla,  KY  40213-3678  . 1 

181  Weather  night,  8160  W  JalTaraon  Blv,  Dallaa,  TX  76211-9670 . 1 

196  Weather  night,  4146  Naval  Air  Rd.,  Pnrt  Huanana,  CA  03041-4001 ...  1 

190  Weather  night,  Wheeler  AFB,  HI  96854-6000  . 1 

200  Weather  night.  6880  Beulah  Rd.,  Sandatnn.  VA  231608100 . 1 

202  Weather  night,  Otla  ANGB,  MA  02642-6028  . 1 

203  Weather  night.  Ft  Indlantawn  GAP,  Annvilla,  PA  17003-6(X)t . 1 

204  Weather  night.  McGuire  AFB,  NJ  08841-6004 . 1 

207  Waathar  night.  3668  N.  Mlehigwi  Rd.,  ShtlbyviUa,  IN  48178-4014  ...  1 

208  Weather  night,  306  AIrpoK  DB,  St  Paul,  MN  88107-4008  . 1 

200  Weather  night,  PO  Box  6218,  Auatln,  TX  787834218 . 1 

210  Weather  night,  OnUrio  ANGB,  CA  91761-7627 . 1 

COMNAVOCEANCOM.  Coda  N31Z,  StannU  Spacn  Ctr.  MS  3S6I04OOO  ...  2 
COMNAVOCBANCOM  (Capl  Brown,  Coda  N332X  Stannla  Spacn  Or. 

MS  30629-6001  . 1 

NAVOCEANO  (Ruaty  RuaaumX  Bldg  8100,  Rm  2030,  Stannla  Space  Cirv 

MS  30622-8001  . 2 

NAVOCEANO,  Code  0230  (Tony  OctolanoX  Stannla  Space  Or;  MS 

30620-6001  . 1 

Maury  Oceanographic  Library,  Naval  Oceanography  Offlee. 

Stannla  Space  Ctr.  MS  30622-8001  . 1 

Naval  Raaaarch  Labomtocy,  Montarey.GA  93043-8006  . . ........I 

Naval  Raaaarch  laibcralory.  Coda  4333,  Waahktgtan.  DC  80378  .........  1 

Naval  natgraduata  School,  Chmn,  Dapt  of  Meteorology,  Coda  83k  Mtetatey,. 

CA  93943-6000  . 

Naval  Eaatam  Ooaanagraphy  Ctr  (Cllm  Sact)cn)i,  U117  McCAdji  BMo, 

Norfolk  NAS,  Nocfblk,  VA  23611-6000 . .1 
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N«»«l  WMtoin  OcBUOfnvhy  Ctr,  Bra  1 13,  AUn:  T«eh  Ubraiy,  PmH  Hwtwr, 

HI9S8aO«)00 . . 

Coramandlnc  Omeer,  Naval  Mar  Oeaanoiraphy  Cmlar,  4301  Sultland  Road, 

FOB  N.  Waahiatton.  DC  30306-6108  . 1 

Naval  Oeaanotnphy  Command  Ctr.COMNAVMAR  Bra  12,  FPO  San  Franoiaeo, 

CA  900306000 . 1 

Commandini  Oflkar,  Naval  Ooeanofraptiy  Command  Ctr,  PSC  819,  Bra  13, 

pro  AS,  0e046-3900 . 1 

NAVOCEANOOMDCT,  Fadaral  BuUdini,  Aahavllle,  NC  28801-2000 . 1 

NAVOCEANCOMDBT,  PaUuent  River  NAS,  MD  20670-5103  . I 

NAVOCBANOOMPAC,  NAS  Norih  laland,  San  Dtofo,  CA  921366130  _ 1 

Naval  Air  Warfare  Cantar-Waapona  DIviaion,  Gaophytkal  Scianeea  Branch, 

Coda  3264,  AUn:  Mr.  Roger  Helvey,  Pbint  Mufu,  CA  930426001  . 1 

WSO,  H  *  HS  Mmtoa  SUtion  Wea,  MCAS  Toetln  CA  02710-6000  . 1 

Araiad  Fonaa  Medical  InUUIgence  Center,  Information  Servicea  Diviiion 

Bldg  1007,  FI  Detriek,  Pradarldi,  MD  217026004  . I 

TECOM,  Dir  for  Tachnotogy,  APO  Met  Team,  Bldg  1134, 

AM8TB-TC-AM  CAB,  Aberdeen  Proving  Ground,  MO  210066001  . I 

Atmoapherie  Sciencea  Laboratory  (SLCAS-AS-I  3  10-20), 

White  Sanda  Miaaile  Range,  NM  880026601 . 1 

TECOM  Atmoa  8d  DIv,  AM8TE-TC-AA  (MacBlain), 

White  Sanda  Miaaila  Range,  NM  880026604 . 1 

White  Sank  Mel  Tbam.  AMSTE-TC-AM  (WSX 

White  Sanda  Miaaile  Range,  NM  880026601 . 1 

Army  Miaaila  Command,  ATTN:  AMSMI-RO-TE-F, 

Radatone  Araanal,  AL  368006260 . 1 

USATBCOM,  ATTN;  AMSTE-TC-AM  (RE)  TCOM  Met  Team, 

Radatene  Araanal,  AL  36808-8062  . I 

USATBCOM,  ATTN;  AMSTE-TC-AM  (AB),  Aberdeen  Proving  Ground, 

MD  210066001 . 1 

Director,  U3.A.-CETEC,  AUn:  GL-AE  (Whitmareh),  PoK  Belvoir, 

VA  220006648 . 1 

Technical  Library,  Dugeray  Proving  Ground,  thigway,  tTT  84022-6000  ....  1 


OPCM,  Attn:  Col  Ihimant,  Suita  900,  0010  Eaeeutive  Blvd. 

Rockville,  MD  20062  . I 

HQ  NATO  Stair  Metaorologleal  Officer  IMSCP8  APO  AE09T34 . 1 

NOAA/MA8C  Library  MC6,  325  Broadway,  Bouldar,OO80a(»6a28 . X 

NOAA  Ubrary-EOC4W6C4,  AUn:  ACQ,  0009  Eaeeutive  Blvd.  RockvIBa. 

MD  20862  . 1 

NOAA/NESDIS  (AUn:  Nancy  Evereon,  E(RA22),  World  Weatbar  Bldg, 

Rm  703,  WaahIngton,  DC  20233  . 1 

NGDC,  NOAA,  MaU  Code  B43C4.  326  Broadway,  Boulder,  GO  803336308 . 1 
NWS  W/OSD,  Bldg  S8M  C-2  Eaat-Weet  Hwy,  SUver  Spring,  MD  20910  ...  I 

NWSTmlnlngCanlar,ei7Hardeaty,Kanaaa  City.  MO  84124  . 1 

NIST  Puba  Production,  Rm  A63e,  Admin  Bldg.  Gailbaraburg,  MO  20899  . .  1 

DTIC-PDAC,  Cameron  Station,  Alaaandria,  VA  223046146  . 2 

AUIASE,Maawell  APB,  AL  361 126564  . 1 

AWSTU  ScoU  APB,  IL  02226-6438  .  36 
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